LABORATOIRE DE PHYSIQUE THEORIQUE ET HAUTES ENERGIES 



THESE DE DOCTORAT DE L'UNIVERSITE 

PARIS 6 



Speciality : PHYSIQUE THEORIQUE 

presentee par 

Julien SIEBERT 



pour obtenir le grade de 
Docteur de l'Universite Paris 6 



Sujet : 

Mecanique statistique 
des gaz autogravitants 



Soutenue le 24 juin 2005 devant le jury compose de : 



M. Hector DE Vega, directeur de these 

M. Claudio Destri, rapporteur 

M me Norma Sanchez , rapporteur 

M. Vladimir Dotsenko, examinateur 



ii 



mes grands-parents 



iv 



Remerciements 



Je suis reconnaissant a Laurent Baulieu de m'avoir accueilli au LPTHE dont il etait 
le directeur et a Hector de Vega d'avoir dirige mon travail de these avec beaucoup de 
competence. Je me felicite que Norma Sanchez, Claudio Destri et Vladimir Dotsenko 
aient accepte d'etre membres de mon jury de these ; nous avons eu des discussions tres 
interessantes et tres fructueuses. J'exprime ma gratitude a Valerie Sabouraud pour l'aide 
administrative qu'elle m'a apportee. Je remercie tous les professeurs qui m'ont fait aime 
la physique. 

Je remercie tous ceux qui m'ont soutenu pendant mon doctorat. Mes pensees vont a 
ma famille, a ma soeur Laurence et surtout a mes parents qui sont si chers a mon coeur. 
Je salue mes amis et particulierememt, parmi eux, Gael Lapeyronnie, Florence Renard, 
Nicolas Gallaud, Roger Gottlieb et Jean-Frangois Richard. Je n'oublie pas de remercier 
mon ami, le docteur en medecine, Christian Lejeune. J'ai egalement des penses amicales 
pour Jonathan Ortiz et Gwendoline Petroffe. 



v 



vi 



Table des matieres 



Introduction 



Hvdrostatiauel 



1.1 



Exemplea 



1.4 Gra nd^ursjjhrokjue'i 



L4J^^je_ 2a£ametre rj 



lA2__L^eq jm^n d'etat 



1.4.3 L'energid 



1.4.4 



1.4.5 



1.5 Cas limited 



_Ch^lejrrs_s_2edfiq uesl 
Compressibilu^sl . . 



11.5.1 p^gyglgg Dement ^ jmu^^tem^era iurel 

11.5.2 Developpement autour de la sphere singulicrd 



1.6 Stabilitd 



1.6.1 Ensemble microcanoniaue et ensemble canoniquel 



L^^^^^bili^^^^ ^hir^^U^rmer 

1.6.3 Instabilites gravitationnelles en astroplivsiqud 



2 Meca jnajjg_statisticju3 

12.1 Ensemb le_jrncjxKaJlQnifli id . . 
l2J_J_^ombj|e_d^jr]jCToet at d 
12.1.2 Grandeurs phvsiaued . 



2.2 Ens ^mbJc_£anoniqjj.d 



l2 : 2J^^bnction_d^_ 2artition 
12.2.2 Grandeurs phvsiqucd 



2.3 Champ movcnl 



l2J3J^_L^enger^ . . . . 

I^^^^^ejj^emble^r iicrocanoniqud 

2. I Calcnls Monte Carlo! 



2.4.1 Algorithme de Metropolid 



2.4.2 R gguitat g dans l'e j^embJe_£anoniguel 



2.4.3 Resultats dans l'ensemble microcanoniquj 



Vll 



3 Svstemes autogravitants comportant plusieurs sortes de particulesl 

3.1 Hvdrostatiaue! 



3.1.1 Eauilibre thermodvnamiauel 



3.1.2 Grandeurs physiques! 



3.1.4 Lois d'echelle 



3.2 



M ccanici ue st at istia uc 



3.2.1 Fonction de partition 

3.2.2 Grandeurs nhvsicmesl 



3-2.3 Champ moven 



4 Svstemes autogravitants en presence de la constante cosmologique! 



4.1 La constante cosmologique 



4.1.1 La constante cosmologique et l'energie noirel 



.4.1.2 Gravitation non relativiste 
4.2 Hvdrosta,tiquel 



4.2.1 Eauilibre hvdrostatiauc 



4.2.2 Eauilibre thermodvnamiouel 

4.2.3 Les parametres r? et fl . . 



4.2.4 Densite de la sphere isotherme 



4.2.5 Stabilite de la sphere isothermel 



4.3 Mccaniauc statistiauc 



4.3.1 Fonction de partition 



4.3.2 Champ movenl 



4.3.3 Calculs Monte Carlo! 



4.4 Limitc R\ ^> 1 



4.5 Discussion; 



Conclusions et perspectives! 



Appendices! 



A Gaz autogravitants polvtropiauesl 



A.l Les transformations polytropiquesl 
A. 2 Lesetoiled 



A. 3 Equilibre hvdrostatiqucl 



B Theorie de Jeans! 



B.l Equation de propagation 
B.2 Relation de dispersion! 



B.3 Instability dans les fluides homogenes isothermesl 
Bibliographie! 



viii 



Introduction 



La formation de structures dans l'univers est gouvernee par la gravita- 
tion (H [2] . Les systemes autogravitants sont des systemes de particules qui 
interagissent entre elles par la gravite ; ils peuvent decrire les distributions 
de galaxies et le milieu interstellaire et notamment expliquer leur structure 
autosimilaire [3J . Les etoiles sont aussi decrites par les systemes autograv- 
itants (avec des effets quantiques) [U, [S]. Nous allons nous interesser aux 
systemes autogravitants thermalises et nous allons presenter la mecanique 
statistique des systemes autogravitants (ne comportant que des particules 
identiques) 0, [7|, |K| pour decrire ensuite des systemes autogravitants avec 
des particules differentes. 

Nous sommes habitues a etudier les systemes de particules en interac- 
tion a courte portee. Ces systemes sont homogenes a l'equilibre thermo- 
dynamique. Leur energie et d'autres grandeurs thermodynamiques comme 
leur entropie, leur energie libre,... sont extensives, c'est a dire qu'elles sont 
proportionnelles au nombre de particules N du systeme. Dans la limite 
thermodynamique usuelle ou le nombre de particules N et le volume V 
tendent vers l'infini et ou le rapport y est fini, les ensembles microcanon- 
ique, canonique et grand-canonique donnent les memes resultats physiques. 
Comme consequence de l'interaction a longue portee, les systemes au- 
togravitants ont des proprietes differentes. Ils ne sont pas homogenes a 
l'equilibre thermodynamique quoique leur energie, leur entropie, leur energie 
libre,... soient extensives (elles sont proportionnelles au nombre de partic- 
ules N du systeme). La limite thermodynamique N — > oo, V — > oo avec 
-T- fini est la limite thermodynamique pertinente pour decrire les systemes 

V3 _ 

autogravitants thermalises car dans cette limite l'energie thermique dont 
la valeur caracteristique est NT et l'energie autogravitationnelle dont la 
valeur caracteristique est Gn ^ 1 N sont du meme ordre (G est ici la constante 

de gravitation, m la masse des particules et T la temperature). Lorsque 



1 



l'energie thermique domine l'energie autogravitationnelle, le systeme se 
comporte comme un gaz parfait. Lorsque l'energie autogravitationnelle 
domine l'energie thermique, le systeme collapse sous l'effet de l'autogravite ; 
c'est d'ailleurs ce qui arrive pour ces systemes dans la limite standard 
N — > oo, V — > oo avec y fini. La limite thermodynamique du gaz auto- 
gravitant N — ► oo, V — > oo avec — r fini est une limite thermodynamique 

diluee ; en effet, la densite moyenne y tend vers comme N~ 2 avec N — > oo. 
Pour les systemes autogravitants les ensembles statistiques ne sont pas 
equivalents ; dans l'ensemble microcanonique, un systeme autogravitant 
peut exister avec une chaleur specifique negative alors que dans l'ensemble 
canonique, il doit avoir une chaleur specifique positive. 

Les systemes autogravitants thermalises peuvent exister sous deux pha- 
ses qui ne peuvent pas coexister, une phase gazeuse et une phase collapsee. 
Dans l'approche du champ moyen, la fonction de partition est evaluee 
pour N — > oo par une integrale fonctionnelle sur la densite ; le poids 
de chaque densite dans l'integrale fonctionnelle est l'exponentielle d'une 
"action" effective proportionnelle a N. On applique 1' approximation de 
point col en cherchant la densite de point col qui rend Taction effective 
extremale. Lorsque Taction effective est minimale, la densite de point col 
domine l'integrale fonctionnelle et l'approche du champ moyen ainsi definie 
decrit exactement la phase gazeuse dans la limite thermodynamique auto- 
gravitante (N — > oo , V — > oo avec -^r fini). La densite de point col est 
solution d'une equation aux derivees partielles du second ordre. En fait, 
la densite de point col correspond a la densite d'un fluide autogravitant a 
Tequilibre hydrostatique obeissant localement a Tequation d'etat des gaz 
parfaits : en chaque point q du systeme, la pression P(q) et la densite de 
masse p m (q) sont liees par la relation 

P(q) = - Pm (q) (1) 

m 

ou T est la temperature constante et m est la masse des particules. Le 
systeme se comporte comme un fluide autogravitant en equilibre hydro- 
statique, les forces gravitationnelles engendrees par l'ensemble du systeme 
autogravitant et les forces de pression se compensant en chaque point ; 
l'approche hydrostatique des systeme autogravitants |H El HOI EH D2i 
fTHl [Hj se deduit done de la mecanique statistique dans l'approche du 
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champ moyen. La mecanique statistique dans l'approche du champ moyen 
determine l'equation d'etat (jlj) alors que dans l'approche hydrostatique, 
l'equation d'etat n'est pas determined, elle doit etre supposee. II est raisonnable 
que le systeme autogravitant obeisse localement a l'equation d'etat des 
gaz parfaits dans cette limite diluee. En symetrie spherique, l'equation 
de la densite devient une equation differentielle du second ordre, appelee 
equation de Lane-Emden. Chacune de ses solutions correspond a une con- 
figuration d'equilibre qui depend des conditions imposees au systeme (tem- 
perature T, rayon Q de la paroi spherique qui enferme le systeme, pression 
P exercee sur la paroi). On peut, par exemple, representer ces configura- 
tions d'equilibre sur le diagramme de phase (fig Jl.ll) ou est tracee la courbe 
f(rj R ) avec 

PV 

et 

7? Cm 2 iV 

Sur ce diagramme, chaque point represente une configuration d'equilibre, 
le point (j) = 0, / = 1) representant la limite des hautes temperatures ou 
l'autogravite (les interactions mutuelles entre particules) est negligeable 
devant l'agitation thermique et ou le systeme se comporte comme un gaz 
parfait (homogene). Pour les autres configurations, / est inferieur a 1 car 
l'autogravite a pour effet d'attirer les particules vers le centre, ce qui fait 
diminuer la pression sur la paroi. Les grandeurs thermodynamiques s'ex- 
priment en fonction de r] R et de /. Les configurations d'equilibre sont iden- 
tiques dans l'ensemble microcanonique, dans l'ensemble canonique et dans 
l'ensemble grand-canonique, mais leur stabilite est differente dans ces trois 
ensembles. Ceci vient du fait que les contraintes imposees au systeme sont 
differentes dans ces ensembles ; l'energie est fixe dans l'ensemble micro- 
canonique mais n'est pas fixe dans l'ensemble canonique. Dans l'ensemble 
canonique (respectivement microcanonique), les configurations d'equilibre 
stable ont leur compressibilite isotherme (respectivement adiabatique) pos- 
itive et sont compris sur le diagramme de phase (fig Jl.ll) entre le point 
(77 = 0, / = 1) et le point ou la compressibilite isotherme (respectivement 
adiabatique) diverge en devenant negative. Les calculs Monte Carlo mon- 
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trent que le champ moyen decrit la phase gazeuse avec une grande precision. 
Dans l'ensemble canonique (respectivement microcanonique) , la transition 
de la phase gazeuse vers la phase collapsee s'opere lorsque la compressibilite 
isotherme (respectivement adiabatique) diverge en devenant negative, en 
conformite avec les previsions du champ moyen. La transition de phase 
se traduit par une grande discontinuite de la pression avec une pression 
negative dans la phase collapsee. 

Le milieu interstellaire est constitue de plusieurs sortes d'atomes (hy- 
drogene,helium,...) et de molecules (dihydrogene,monoxyde de carbone,...) 
et les distributions de galaxies sont constitutes de galaxies de masses 
differentes. Nous avons done etudie la mecanique statistique des systemes 
autogravitants comportant n sortes de particules (iV~i particules de masse 
?rzi, N2 particules de masse m2,..., N n particules de masse m^) jT5]. La lim- 
ite thermodynamique pertinente est telle que les nombres de particules Ni 
et le volume V tendent vers l'infini et les rapports -^f- sont finis. Nous 
avons developpe l'approche du champ moyen de la mecanique statistique 
qui decrit exactement la phase gazeuse dans cette limite thermodynamique. 
La mecanique statistique dans l'approche du champ moyen montre que le 
systeme se comporte comme un melange de gaz autogravitants differents 
obeisssant chacun localement a l'equation des gaz parfaits : en chaque point 
q, la pression partielle Pi et la densite de masse partielle p m% des particules 
de masse m 2 sont liees par la relation 

Pi(q) = ^-PmM ■ (4) 

Chaque gaz autogravitant est en equilibre hydrostatique (en chaque point, 
les forces gravitationnelles exercees sur les particules de masse m z par 
la totalite du systeme et les forces de pression exercees sur les partic- 
ules de masse rrii par les particules de masse rrii voisines se compensent). 
Les densites partielles sont solutions d'un systeme de n equations aux 
derivees partielles du second ordre qui se reduit a une seule equation aux 
derivees partielles du second ordre et qui en symetrie spherique devient 
une equation differentielle du second ordre. Chacune de ses solutions cor- 
respond a une configuration d'equilibre qui depend des conditions imposees 
au systeme (temperature T, rayon Q de la paroi spherique qui enferme le 
systeme, pression partielle Pi exercee sur la paroi par les particules de 
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masse mi). Ces configurations d'equilibre sont representees par les points 
du diagramme de phase constitue par les courbes fi(r)i , N2, N n ) 
fi(N u Ni-i, rf , Ni+i, N n ) , f n (Ni, ^ ) ou 

et 

^ = ^i. ( 6 ) 

L'energie, l'entropie, l'energie libre, les chaleurs specifiques a pression cons- 
tante et a volume constant ainsi que les compressibilites isothermes et adi- 
abatiques sont calculees en fonction de fi, f n et de 7/f , rj^. Nous avons 
trouve que ces grandeurs thermodynaniques sont extensives ; dans la limite 
Ni — > 00, A^2 — >• cxd, — >• 00, chaque grandeur s'exprime comme la 
somme d'un terme proportionnel a Aq, d'un terme proportionnel a A^2, ... 
et d'un terme proportionnel a N n . Par l'etude des densites partielles, nous 
avons montre que les particules les plus lourdes subissant plus intensement 
les effets attractifs de l'autogravite sont plus nombreuses pres du centre de 
la sphere tandis que les particules les plus legeres sont plus nombreuses pres 
de la paroi. Nous avons analyse la stabilite de ces systemes dans l'ensemble 
canonique et dans l'ensemble microcanonique. Nous avons montre que les 
systemes autogravitants comportant deux sortes de particules obeissent a 
des lois d'echelle liant leur masse M et leur taille q 

M{q) - q d 

ou la dimension fractale d est inferieure ou egale a 3. La dimension fractale 
d est en generale dependante de la composition du melange (c'est a dire 
de j^-). Cependant aux points critiques du systeme ou la chaleur specifique 
a volume constant diverge, la dimension fractale est independante de la 
composition du melange et vaut 1.6.... Ceci manifeste 1' "universalite" des 
proprietes du systeme a l'approche du regime critique. 

Les recentes observations astrophysiques ont montre que l'univers est 
rempli de ce qui est appele l'energie noire et qu'on modelise par la cons- 
tante cosmologique des equations de la relativite generale [TBI El HBj • Elle 
agit comme une densite d'energie uniforme ayant une action gravitation- 
nelle repulsive sur la matiere. II est done important d'etudier son influence 
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sur les systemes autogravitants. Nous avons etudie la mecanique statistique 
des systemes autogravitants (une seule sorte de particules) en presence de 
la constante cosmologique [191 ED|- Nous avons trouve que la limite thermo- 
dynamique pertinente des systemes autogravitants en presence de la cons- 
tante cosmologique est la limite ou le nombre de particules N et le volume 
V tendent vers l'infini et ou la constante cosmologique A tend vers avec 
-^r et AV^ finis. L'energie de la constante cosmologique dans le systeme 
a pour valeur caracteristique AV, tandis que l'energie thermique a pour 
valeur caracteristique NT et l'energie autogravitationnelle a pour valeur 
caracteristique Gn ^i I ■ Pour que l'energie de la constante cosmologique soit 
de meme ordre que l'energie thermique il faut que AV soit de l'ordre de N 
lorsque N — > oo ; or le rapport -^r est fini lorsque N — > oo et V — > oo pour 
que l'energie thermique et l'energie autogravitationnelle soient de meme 

2 

ordre. La quantite AV s doit done etre finie lorsque A — ► et V — > oo pour 
que l'energie de la constante cosmologique, l'energie autogravitationnelle 
et l'energie thermique soient de meme ordre. Dans cette limite thermody- 
namique diluee, l'approche du champ moyen decrit exactement la phase 
gazeuse et montre que le systeme se comporte comme un gaz en equilibre 
hydrostatique obeissant localement a l'equation des gaz parfaits (CD). En 
chaque point, les forces gravitationnelles attractives exercees par l'ensem- 
ble du systeme autogravitant, les forces gravitationnelles exercees par la 
constante cosmologique et les forces de pression se compensent. La densite 
de masse est solution d'une equation aux derivees partielles du second ordre 
qui en symetrie spherique devient une equation differentielle du second or- 
dre. Chaque solution correspond a une configuration d'equilibre dependant 
des conditions imposees au systeme (temperature T, rayon Q de la paroi 
spherique qui enferme le systeme, pression P exercee sur la paroi par les 
particules, constante cosmologique A). Les configurations d'equilibre sont 
representees sur le diagramme de phase (fig J4.ll ) ou sont tracees les courbes 
f(rj R , Ra)- Les parametres r] R et / sont definis par les equations (j2j) et (jHj) ; 
Ra est le rapport entre l'energie de la constante cosmologique et la masse 
de la matiere 

* - ■ 

La constante cosmologique exerce une action gravitationnelle repulsive sur 
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la matiere qui s'oppose a Taction attractive de l'autogravite. En absence 
de la constante cosmologique (R\ = 0), la densite comme fonction de 
la distance q par rapport au centre de la sphere (0 < q < Q) est une 
fonction decroissante. En presence de la constante cosmologique, la densite 
peut etre soit decroissante, soit uniforme, soit croissante. Lorsque R\ < 1, 
les effets attractifs de l'autogravite l'emportent sur les effets repulsifs de 
la constante cosmologique ; la densite est une fonction decroissante de q. 
Lorsque R\ = 1, les effets attractifs de l'autogravite et les effets repulsifs 
de la constante cosmologique se compensent exactement ; la densite est une 
fonction uniforme de q et le systeme se comporte comme un gaz parfait (la 
pression P est uniforme et verifie la loi des gaz parfaits jtL = 1). Lorsque 
R\ > 1, les effets repulsifs de la constante cosmologique l'emportent sur 
les effets attractifs de l'autogravite ; la densite est une fonction croissante 
de q. Nous avons calcule les grandeurs thermodynamiques ; nous avons 
trouve qu'elles sont extensives (proportionnelles au nombre de particules 
N). Nous avons determine la zone de stabilite de ce systeme suivant R\. 
Nous avons fait des calculs Monte Carlo pour ce systeme [20] • Ces calculs 
Monte Carlo dans l'ensemble canonique montrent que le champ moyen 
decrit tres bien la phase gazeuse et que la transition de la phase gazeuse 
vers la phase collapsee s'opere lorsque la compressibilite isotherme diverge 
en devenant negative. Nous avons egalement etudie le cas limite Ra ^> 
1 ou l'autogravite est negligeable devant la constante cosmologique. Les 
particules sont soumises a des forces harmoniques dont la pulsation au 
carre est negative. Dans ce cas limite, le systeme est exactement soluble ; 
nous avons calcule les grandeurs thermodynamiques et nous avons trouve 
que toutes les configurations d'equilibre sont stables. 

Dans les deux premiers chapitres, nous presentons les systemes auto- 
gravitants thermalises constitues par une seule sorte de particules. Dans 
le premier chapitre est exposee l'hydrostatique des systemes autogravi- 
tants. Nous etablissons l'equation verifiee par la densite, nous calculons les 
grandeurs thermodynamiques et nous explorons la stabilite du systeme. 
Dans le deuxieme chapitre, la mecanique statistique des systemes auto- 
gravitants dans l'ensemble microcanonique et dans l'ensemble canonique 
est presentee. Nous montrons que dans l'approche du champ moyen, l'hy- 
drostatique est deduite et que l'equation d'etat est derivee ; il s'agit de 
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1' equation locale des gaz parfaits inhomogenes. Nous exposons enfin les 
calculs Monte Carlo. 

Dans le troisieme chapitre, nous presentons les resultats obtenus sur 
les systemes autogravitants constitues de plusieurs sortes de particules. 
Dans l'approche du champ moyen, nous derivons les equations d'etat, nous 
etablissons les equations des densites partielles, calculons les grandeurs 
thermodynamiques, faisons 1' analyse de la stabilite du systeme et etudions 
ses lois d'echelle. 

Dans le quatrieme chapitre, nous presentons les resultats obtenus sur les 
systemes autogravitants en presence de la constante cosmologique. Dans 
l'approche du champ moyen, nous derivons l'equation d'etat, etablissons 
1' equation de la densite et etudions la stabilite des configurations d'equilibre 
dans 1' ensemble canonique. Nous exposons le cas limite R\ ^> 1. Nous 
presentons les calculs Monte Carlo. 
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Chapitre 1 
Hydrostatique 



Un systeme autogravitant est un systeme de particules qui interagissent 
entre elles par la gravitation. On se limite ici a la gravitation newtonienne. 
II existe sous deux phases qui ne peuvent pas coexister ensemble, une phase 
gazeuse et une phase collapsee. L'outil adequat pour etudier les proprietes 
de ce systeme a l'equilibre thermodynamique est la mecanique statistique 
(chapitre 2). L'approche du champ moyen dans laquelle le systeme se com- 
porte comme un gaz autogravitant isotherme en equilibre hydrostatique et 
obeit localement a l'equation d'etat des gaz parfaits decrit exactement la 
phase gazeuse dans la limite thermodynamique ou le nombre de particules 
N et le volume V sont tels que N — > oo, V — > oo et -\ est fini. Dans la 

limite thermodynamique standard ou N — > oo, V — > oo et y est fini, le 
systeme est dans sa phase collapsee. Nous allons etudier dans ce chapitre 
l'hydrostatique des gaz autogravitants parce qu'elle decrit exactement la 
phase gazeuse dans la limite thermodynamique. Nous allons voir que Ton 
peut appliquer le gaz autogravitant a l'equilibre thermodynamique a au 
moins deux types d'objets astrophysiques : le milieu interstellaire et les 
distributions de galaxies. 

1.1 Exemples 

Le gaz autogravitant a l'equilibre thermodynamique decrit tous les ob- 
jets astrophysiques thermalises dans lesquels les interactions gravitation- 
nelles jouent un role preponderant. Le milieu interstellaire est constitue 
de nuages de gaz et de poussieres situes dans le plan des galaxies. C'est 
dans ces nuages que se forment les etoiles par effondrement gravitationnel. 
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Le milieu interstellaire n'est pas homogene, les nuages le formant ont des 
tailles diverses, leur taille etant inversement proportionnelle a leur densite. 
Les nuages les moins denses contiennent de l'hydrogene atomique HI, les 
plus denses contiennent de l'hydrogene moleculaire H^ et des molecules 
formees a partir d'elements plus lourds comme le monoxyde de carbone 
CO. Le milieu interstellaire a une structure tres interessante. En effet, les 
observations des raies moleculaires [2H [22] ont permis d'obtenir des in- 
formations sur la masse et la dynamique de ces nuages. II a ete mis en 
evidence, pour des regions ayant une longueur / comprise entre 1CT 2 par- 
sees et 100 parsecs que la dispersion interne des vitesses Av et la masse M 
obeissent a des lois de puissance 

M~l dH , Av~l dv (1.1) 
avec des dimensions fractales d# et d v verifiant 

1-4 < d H < 2 , 0.3 < d v < 0.6 . 

De telles relations montrent que le milieu interstellaire a une strucure au- 
tosimilaire se repetant a toutes les echelles [23J . Cette strucure est hierarchi- 
que, les grands nuages sont fragmentes en de plus petits nuages condenses, 
ceux-ci etant fragmentes en de plus petits nuages encore plus condenses et 
ainsi de suite sur au moins 5 a 10 ordres de grandeurs. La limite superieure 
de cette clusterisation est de 100 parsecs, ce qui correspond a un million de 
masses solaires. Des nuages de taille superieure ne peuvent pas exister car 
ils seraient detruits par les forces de marees galactiques. La limite inferieure 
est limitee par la resolution des telescopes, il semble qu'elle puisse etre re- 
poussee jusqu'a 10~ 4 parsecs, ce qui correspond a une masse de l'ordre de 
celle de Jupiter. L'autosimilarite du milieu interstellaire sur une echelle de 
tailles aussi larges a suscite de nombreuses tentatives d'explication [22j. 
La physique du milieu interstellaire est complexe (formation d'etoiles, ex- 
plosion de super novae, vent stellaire...). On peut cependant poser comme 
hypothese que l'essentiel de la physique du milieu interstellaire vient des 
interactions gravitationnelles entre les particules qui le composent. On peut 
supposer que le milieu interstellaire est thermalise sur de larges echelles. 
Le systeme autogravitant isotherme est done utile et plein de sens pour 
la description du milieu interstellaire |3j. Or, il a ete montre que les gaz 
autogravitants a l'equilibre thermodynamique obeissent a des lois d'echelle 
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sur la masse et la vitesse analogues a celles de la relation ( 11.11 ) [7j . Les gaz 
autogravitants a l'equilibre thermodynamique sont done aptes a decrire et 
expliquer la structure autosimilaire du milieu interstellaire. Les gaz auto- 
gravitants s'appliquent aussi aux distributions de galaxies, chaque galaxie 
etant considered comme un point materiel interagissant avec les autres 
galaxies par l'intermediaire de la gravite ; cependant, les galaxies ne sont 
pas thermalisees. Comme le milieu interstellaire, les galaxies s'organisent 
en structure hierarchique du groupe de galaxies qui est la structure la plus 
petite, en passant par l'amas de galaxies jusqu'au superamas de galax- 
ies qui est la structure la plus grande. La structure des distributions de 
galaxies est autosimilaire et leur masse et taille sont reliees par une loi de 
puissance analogue a celle du milieu interstellaire ( 11.11) avec une dimen- 
sion dn ~ 1.7 et ce pour une echelle allant jusqu'a 10 9 annees lumieres 
[24|. Au dela, la dimension fractale tend vers 3 car l'expansion de l'univers 
l'emporte sur l'autogravite et l'univers devient homogene. Le fait que le 
milieu interstellaire et les distributions de galaxies aient le meme type de 
structure autosimilaire, avec les memes lois d'echelle, plaide en faveur de la 
gravite comme cause de ce phenomene universel. Ceci motive tres fortement 
l'etude des gaz autogravitants. Nous allons presenter les gaz autogravitants 
en equilibre hydrostatique. 

1.2 Equilibre hydrostatique 

Dans l'approche hydrostatique, le gaz autogravitant est divise en ele- 
ments de fluide, chaque element de fluide subissant deux forces qui se com- 
pensent : 

-les forces gravitationnelles exercees par la totalite du gaz autogravitant 
-les forces de pression exercees par les elements de fluide qui lui sont 
voisins. 

Le gaz autogravitant est alors en equilibre hydrostatique ; l'equilibre des 
forces de pression et des forces gravitationnelles s'appliquant sur un element 
de volume unite centre autour du point q s'ecrit 



ou P est la pression du gaz au point q, p m est la densite de masse au point 



V*P + p m (q) 9(q) = 
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q et 

M = " G [ d 3 q'p m (q') (1.3) 

J \q-q \ 

est le champ gravitationnel engendre au point q par le gaz autogravitant. 
En utilisant l'equation de Poisson du champ gravitationnel 



= -4rcG Pm (q) , (1.4) 

on obtient la relation suivante qui est la condition d'equilibre hydrostatique 
d'un gaz autogravitant 

V«f ( — V q P) = -4ttG p m {q) . (1.5) 

Nous allons souvent considerer les gaz autogravitants en symetrie spheri- 
que que nous allons appeler spheres autogravitantes. A cause de la symetrie 
spherique, toutes les grandeurs physiques ne dependent que de la distance 
q par rapport au centre de la sphere. Soit M(q) la masse a l'interieur de la 
sphere de rayon q : 

M(q) = du 4tt u 2 p m {u) . (1.6) 
Jo 

En symetrie spherique l'integration du champ ( 11.31) donne 

G M{q) p m {q) 



g(q) = 



Q 2 



En utilisant l'equation ( 11.21) . on a 

dP G M(q) p m (q) 



dq q 2 

On remarque d'apres cette relation qu'une sphere autogravitante en equi- 
libre hydrostatique a une pression qui decroit en s'eloignant du centre 
(if < ^ es f° rces de pression sont done orientees vers l'exterieur, elles 
s'opposent aux forces gravitationnelles qui tendent a contracter la sphere 
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autogravitante vers son centre. La condition d'equilibre hydrostatique ( 11.51) 
devient en symetrie spherique 




= -47rGp m (q). (1.7) 

Cette relation exprime la condition d'equilibre hydrostatique d'une sphere 
autogravitante. L'equation ( 11.71 ) est la limite newtonienne de la condition 
d'equilibre hydrostatique d'un fluide autogravitant en relativite generale 
(voir par exemple |26j). 

Nous allons presenter maintenant les gaz autogravitants a l'equilibre 
thermodynamique ; ce sont les gaz autogravitants isothermes en equilibre 
hydrostatique. 



1.3 Equilibre thermodynamique 

La condition d'equilibre hydrostatique ( 11.51 ) ne sufht pas pour decrire 
un gaz autogravitant. Pour determiner la densite de masse et la pres- 
sion, il faut l'equation d'etat reliant ces deux grandeurs. Nous allons nous 
interesser aux gaz autogravitants a l'equilibre thermodynamique (dans l'an- 
nexe I, on presente un autre type de gaz autogravitant, le gaz autogravitant 
polytropique). En hydrostatique, il faut postuler l'equation d'etat tandis 
que dans l'approche de mecanique statistique du gaz autogravitant, son 
equation d'etat est derivee (chapitre 2). La mecanique statistique, dans 
l'approche du champ moyen, montre que le gaz autogravitant verifie la 
condition d'equilibre hydrostatique ( 11.51 ). dans la limite thermodynamique 
autogravitante ou le nombre de particules N et le volume V tendent vers 
l'infini avec pr fmi. Dans cette limite thermodynamique diluee la den- 
site moyenne y se comporte comme N~ 2 lorsque N — ► oo. L'approche 
de mecanique statistique permet de deriver l'equation d'etat du systeme. 
Elle montre que la pression et la densite de masse du gaz autogravitant 
obeissent, en chaque point q, a l'equation d'etat des gaz parfaits (QJ 

T 

p (q) = —Pm(q) 
m 

ou m est la masse de chacune des particules du gaz qui sont supposees ici 
identiques et ou T est la temperature du gaz isotherme. A partir de cette 
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equation et de l'equation ( 11.51) . on obtient l'equation de la densite du gaz 
autogravitant a l'equilibre thermodynamique 

V|(lnp TO ) + -^— Pm (q) = 0. (1.8) 
En posant pour la densite 




p m = p e x (1-9) 

ou p est une constante et en introduisant le rayon vecteur sans dimension 
A defini par 



(1.10) 

\j 4:71 Lt lib fj 

on trouve l'equation suivante 

V|x + e x(X) = . (1.11) 

Nous verrons dans le chapitre 2 que cette equation correspond au point 
col de la fonction de partition dans l'approche du champ moyen. Dans le 
cas de la symetrie spherique, on appelle le gaz autogravitant en equilibre 
thermodynamique sphere isotherme |H[Sj. L'equation ( 11.111) devient 

hU x2 S) +eX -°- (112) 

Cette equation correspond a l'equation de Lane-Emden isotherme dans 
l'approche hydrostatique jij. On peut poser que p est la densite au centre 
et en deduire la premiere condition initiale 

X (A = 0) = 0. (1.13) 

Pour que l'equation ( 11.121 ) soit reguliere a l'origine, on impose la deuxieme 
condition initiale 

g(A = 0) = 0. (1.14) 

L'equation ( 11.111) est covariante par la transformation d'echelle suiv- 
ante : si x(^) es ^ solution de cette equation alors, etant donne une constante 
C, la fonction 
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X*(X) = X (CX) + 21nC 



est aussi solution de cette equation. Dans le cas de la symetrie spherique, 
cette propriete de covariance s'enonce de cette fagon : si x(A) est solution 
de 1' equation de Lane-Emden isotherme ( 11.121) alors la fonction 



X*(X) = x(CX) + 21nC (1.15) 

est aussi solution de l'equation de Lane-Emden isotherme. Grace a cette 
propriete, on peut deduire toute une famille de solutions a partir d'une so- 
lution. Par exemple, les solutions regulieres (x'(0) = 0) sont engendrees par 
la solution verifiant x(0) = 0. Appliquer cette transformation a la solution 
dont les conditions aux limites sont ( 11.131 ) et ( 11.141) revient a redefinir la 
constante p dans l'equation ( 11.91 ). On va maintenant montrer que les solu- 
tions telles que x(0) est fini verifient necessairement X (0) = jSj- D'apres 
(ITT2D on a 



d 2 (A X ) 
d 2 A 



= -Xe xW 



done 



(%) =U^) = 4 lim [A e*»l . 
VdA; A=0 2 V d^A J x=0 2 A^oL J 

Si x(0) est fini alors on a bien X (0) = 0. 

Dans la section suivante, nous allons exposer le calcul de quelques grandeurs 
physiques de la sphere isotherme qui est le gaz autogravitant a l'equilibre 
thermodynamique en symetrie spherique, a partir des solutions de l'equation 
de Lane-Emden isotherme ( 11.121 ). 



1.4 Grandeurs physiques 

Nous allons presenter le calcul de quelques grandeurs physiques de la 
sphere isotherme. 
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0.2 1 1 1 1 1 1 J2 1 

0.5 1 1.5 2 2.5 ^ 3 

FIG. 1.1 - La courbe f(r] R ) oil / = ^ (voir equation (fTT^jl ) et r/ K = ^# (voir 
equation ()1.17|) ). Le point (rj R = 0, / = 1) correspond a la limite des basses densites et des 
hautes temperatures oil les effets de la gravitation sont negligeables et oil le systeme est 
homogene. Le parametre r] R atteint sa valeur maximale au point (r) R = 2.51..., / = |) ; 
il s'agit du point de temperature minimale oil la chaleur specifique a volume constant 
diverge en devenant negative. Le point (i] R = 2, / = |) correspond a la limite des hautes 
densites oil le systeme est infiniment dense au centre de la sphere. La courbe a la forme 
d'une spirale qui s'enroule autour du point (r] R — 2, f — |). Les configurations comprises 
entre le point (rj R = 0, / = 1) et le point (r/ R = 2.43...,/ = 0.40...) sont stables dans 
l'ensemble canonique. Le point (rj R = 2.43..., / = 0.40...) represents par le symbole +can 
est le point d'instabilite dans l'ensemble canonique. Les configurations comprises entre le 
point (r] R = 0, / = 1) et le point (r) R = 2.14..., / = 0.26...) sont stables dans l'ensemble 
microcanonique. Le point (j] R = 2.14..., / = 0.26...) represents par le symbole +mic est 
le point d'instabilite dans l'ensemble microcanonique. 
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1.4.1 Le parametre rj 

On suppose que le systeme qui est compose de N particules de masse 
m et dont la masse totale est M = mN est contenu dans un volume V. In- 
troduisons le parametre sans dimension 77 = G 7 — = Gm \ — ■ Le parametre 
rj est le rapport entre deux energies caracteristiques d'une particule en 
interaction avec le systeme, G ™i M qui est de l'ordre de son energie gravi- 
tationnelle et T qui est de l'ordre de son energie cinetique. Lorsque 77 tend 
vers 0, son energie cinetique l'emporte largement sur son energie gravita- 
tionnelle et le systeme se comporte comme un gaz parfait ; c'est le cas des 
systemes terrestres usuels. Lorsque 77 est de l'ordre de 1, l'energie gravita- 
tionnelle de la particule est du meme ordre que son energie cinetique et 
la gravite joue un role important dans la physique du systeme. Lorsque 
rj — ► oo, l'energie gravitationnelle de la particule l'emporte largement sur 
son energie cinetique et le systeme s'effondre immediatement sur lui-meme. 
Voila pourquoi la limite thermodynamique pertinente des systemes auto- 
gravitants est TV — > oo, V — > oo avec fini. Pour cette limite, rt est 

V3 

fini et l'autogravite joue un role important. Par contre, la limite thermo- 
dynamique standard N — > oo, V — > oo avec y fini, correspond pour les 
systemes autogravitants a r] — > oo ou le systeme s'effondre sur lui-meme. 
II a ete estime que 77 ~ 1 pour le milieu interstellaire 0, E| . Ceci confirme 
que l'autogravite joue un role tres important dans la physique du milieu 
interstellaire. 

Calculons 77 a partir des solutions de l'equation de Lane-Emden ( 11.121) 
dans le cas de la sphere isotherme. La masse contenue a l'interieur du 
volume V est 



M = mN = / d 3 q ' p m (q 




En utilisant les equations ( 11.81 ) et le theoreme de Green-Ostrogradski sur 
le volume V, on trouve que la masse est proportionnelle a une integrale sur 
la surface de la paroi entourant le volume V 

M = -—L—<fdSV? lnp m . 
4ttG m J 

Dans le cas de la symetrie spherique ou le volume est une sphere de rayon 
Q, le parametre rj R = G q T n = r\ (^") 3 defini dans I'introduction par la 
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relation (J3J) vaut 



n* = -Q(±Q*pj)^. (Lie) 

D'apres la relation ( 11.91 ) et en introduidant le rayon reduit A = v /4 ZL Gmp £ g 
( eq. ( 11.1 Op ), on obtient 

^ = —qjT- = "Ax (A) • (1.17) 
Nous allons maintenant determiner l'equation d'etat de la sphere isotherme. 

1.4.2 L'equation d'etat 

Les equations (JTJ), ( 11.91 ). ( 11.101) et ( 11.171) permettent de deduire la valeur 
du parametre / defini dans l'introduction par la relation (J2J) 

PV 1 A e x(A) 

/= iVT = "37W (1 ' 18) 

ou P est la pression du gaz sur la paroi et V = 4 ^ est le volume du 
systeme. Ce parametre determine l'equation d'etat du gaz sur la paroi. 
Montrons que / comme fonction de 77 obeit a une equation differentielle 
du premier ordre. D'apres les equations ( II . 12H . ( I1.17p et ( 11.181 ). on a 

/ dA _ A ; ' ' 7j B dA _ A j ' 

On en deduit que / comme fonction du parametre 77 obeit a l'equation 
differentielle du premier ordre suivante 

rj R df 3f+r] R -3 , 

T^ = - 3/-1 ; (L19) 

L'equation differentielle ( 11.191) correspond a une equation d'Abel de pre- 
mier type. La condition initiale de cette equation differentielle est f(rj = 
0) = 1, ce qui correspond a la limite des basses densites et des hautes 
temperatures ou la sphere isotherme tend a se comporter comme un gaz 
parfait homogene. Grace a sa propriete de covariance ( 11.151) . l'equation de 
Lane-Emden isotherme qui est une equation differentielle du second ordre 
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a ete reduite en l'equation du premier ordre ( 11.191) d'inconnue f(r] R ). Ce 
resultat a ete obtenu car la grandeur rj R = —Xx {X) definie par la relation 
( 11.171 ) et la grandeur / = — | definie par la relation ( 11.181) sont 

invariantes par la transformation d'echelle ( 11.151) . 

La courbe f(f] R ) (fig JI.ll) constitue le diagramme de phase de la sphere 
isotherme. Chaque point du diagramme represente une configuration d'equi- 
libre, solution de l'equation de Lane-Emden isotherme (eq.( ll.l2f )). Le point 
(rj R = 0,/ = l) correspond a la limite des basses densites et des hautes tem- 
peratures ou les effets de la gravitation sont negligeables et ou le systeme 
est homogene. Le parametre r] R atteint sa valeur maximale au point (77 = 
2.51...,/ = -|) ; il s'agit du point de temperature minimale ou la chaleur 
specifique a volume constant diverge en devenant negative. Le point (77 = 
2,/ = correspond a la limite des hautes densites ou le systeme est in- 
finiment dense au centre de la sphere. Les configurations comprises entre 
le point (rj R = 0,/ = 1) et le point (rj R = 2.43...,/ = 0.40...) sont stables 
dans l'ensemble canonique. Le point (77 = 2.43...,/ = 0.40...) represente 
par le symbole +can est le point d'instabilite dans l'ensemble canonique. 
Les configurations comprises entre le point (77 = 0, / = 1) et le point 
(rj R = 2.14...,/ = 0.26...) sont stables dans l'ensemble micrcanonique. Le 
point (77^ = 2.14..., / = 0.26...) represente par le symbole +mic est le point 
d'instabilite dans l'ensemble microcanonique. 

Determinons la pression a l'interieur de la sphere. Plagons nous a un 
rayon q inferieur au rayon Q de la paroi spherique. D'apres l'equation 
( 11.101 ). a ce rayon q correspond un rayon reduit X q = ^ A. On en deduit 
d'apres les equations (JTJ) et ( 11.91) que la densite de masse et la pression au 
rayon q sont telles que 

m P(q) _ p m {q) =eX ( # A)_ (L2Q) 
Tp p 

Introduisons le contraste qui est le rapport de la pression au centre de la 
sphere et de la pression sur la paroi de la sphere 

C = = e -x(A). n 21) 



P(Q) 

La fonction %(A) solution de l'equation de Lane-Emden isotherme ( 11.121) est 
negative et decroissante. On en deduit que le contraste est une grandeur 
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superieure a 1 et croissante avec A. Chaque configuration d'equilibre de 
la sphere isotherme, solution de l'equation de Lane-Emden isotherme est 
determinee par une et une seule valeur de A, elle est done determined par 
une et une seule valeur du contraste. 

Nous allons maintenant calculer l'energie. 



1.4.3 L'energie 

Pour calculer l'energie, nous allons appliquer le theoreme du viriel [27J. 
Considerons un corps thermalise dont l'energie potentielle d'interaction Ep 
est une fonction homogene d'ordre n des coordonnees des particules et qui 
a comme energie cinetique E c , comme energie totale E = E c -\-Ep et comme 
pression de paroi P. On a la relation suivante 



2 E c -nEp = (n + 2) E c - nE= 3 PV . 

Pour un gaz autogravitant, l'energie potentielle d'interaction est d'ordre 
n — — 1, on a done 



2 E c + E P = E c + E = 3 PV . (1.22) 

En supposant que le gaz est monoatomique, on a E c = | N T. D' apres 
l'equation ( 11.181) . on obtient le resultat suivant pour l'energie 

E 3 (f-l) . (1.23) 



NT V 2, 

Introduisons le parametre e R = G ^F N 2 qui est le quotient de son energie 
E et de 1' energie Gm q N qui est de l'ordre de son energie gravitationnelle. 
D 'apres les equations ( 11.171) et ( 11.231) . on a 

r _ Q E _ o f~\ (1 24 n 

Alors que r) est le parametre pertinent dans l'ensemble canonique, e R est 
le parametre pertinent dans l'ensemble microcanonique. 

Calculons maintenant la chaleur specifique a volume constant et la 
chaleur specifique a pression constante. 
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Fig. 1.2 - La chaleur specifique a volume constant c v , la chaleur specifique a pres- 
sion constante c p , la compressibilite isotherme kt et la compressibilite adiabatique Ks en 
fonction de In A defini par l'equation (|1.1U|) . La chaleur specifique a volume constant est 
positive pour < A < 8.99, diverge pour A = 8.99, est negative pour 8.99 < A < 34.2 et 
s'annule en redevenant positive pour A = 34.2. La chaleur specifique a pression constante 
est positive pour < A < 6.5, diverge pour A = 6.5, est negative pour 6.5 < A < 25.8 et 
s'annule en redevenant positive pour A = 25.8. La compressibilite isotherme est positive 
pour < A < 6.5 et diverge en devenant negative pour A = 6.5, c'est a dire pour la meme 
valeur oil la chaleur specifique a pression constante diverge. La compressibilite adiaba- 
tique est positive pour < A < 25.8 et diverge en devenant negative pour A = 25.8, c'est 
a dire pour la meme valeur oil la chaleur specifique a pression constante s'annule. 
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1.4.4 Chaleurs specif! ques 



Effectuons d'abord le calcul de la chaleur specifique a volume constant 
c v = jj iw)v ^ n u t n i san t les equations ( 11.171 ) et ( 11.231) . on trouve 



f(rj R )-r} R f'(v R )-\ 



c v — 3 

En utilisant l'equation ( 11.191 ). on obtient 



Cv = V-\ + ^ R + J]—{- (1-25) 

Grace aux equations ( 11.171 ) et ( 11.181) . on trace c v (fig J1.2l) en fonction du 
rayon reduit A ( 11.101) . La chaleur specifique a volume constant est positive 
pour < A < 8.99, diverge pour A = 8.99, est negative pour 8.99 < A < 
34.2 et s'annule en redevenant positive pour A = 34.2. 

Effectuons maintenant le calcul de la chaleur specifique a pression cons- 
tante c p = jj (dr)p - 
En utilisant la formule [27j 

r - r - — )Mlv 
p v N fdP\ 

\dv)' 



T 

et les equations ( 11.171) et ( 11.181) . on a 

C > ~ -2 + 41(11 > /(,«) + If fm - 
En utilisant l'equation ( 11.191 ). on obtient 

3 24( 1; fi - 2)/ 

C " = "2 + 12 / + 6/-^ ■ (L26) 

Grace aux equations ( 11.171 ) et ( 11.181) . on trace c p (fig J1.2|) en fonction du 
rayon reduit A ( 11.101) . La chaleur specifique a pression constante est positive 
pour < A < 6.5, diverge pour A = 6.5, est negative pour 6.5 < A < 25.8 
et s'annule en redevenant positive pour A = 25.8. 

Nous allons enfin calculer la compressibilite isotherme et la compress- 
ibilite adiabatique. 
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1.4.5 Compressibilites 



La compressibilite isotherme et la compressibilite adiabatique ont pour 
expression respective 





T ¥ / S 

Elles expriment la variation de la pression par rapport au volume respec- 
tivement a temperature constante et entropie constante. 

Calculons la compressibilite isotherme pour les configurations d'equili- 
bre de la sphere isotherme. En utilisant les equations ( 11.171) et ( 11.181) . on 
trouve que la compressibilite isotherme est egale a 



NT 



1 



V ' f(v R ) + v if'(v R ) 

En utilisant l'equation ( 11.191 ). on a 



NT 



V 



3 

V 



1 + 



R 



6/ -T) 



R 



1.27) 



Grace aux equations ( 11.171 ) et ( 11.181 ) on trace kt (fig J1.2l ) en fonction du 
rayon reduit A ( 11.101) . La compressibilite isotherme est positive pour < 
A < 6.5 et diverge en devenant negative pour A = 6.5, c'est a dire pour la 
meme valeur ou la chaleur specifique a pression constante diverge. 
La compressibilite adiabatique verifie la relation [2~Tj 

c v 

k s = — kt 



c 



p 



En utilisant les equations ( 11.251) . ( 11.261) et ( 11.271 ). on en deduit 
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3 12/ 2 + (2r] R - 11)/+ 1 



V ~ KS f 48 p + (8^ - 38)/ + r) R ' (L28) 

Grace aux equations ( 11.171) et ( 11.181 ). on trace ^5- (fig J1.2l) en fonction du 
rayon reduit A ( 11.101 ). La compressibilite adiabatique est positive pour < 
A < 25.8 et diverge en devenant negative pour A = 25.8, c'est a dire pour 
la meme valeur ou la chaleur specifique a pression constante s'annule. 

Nous allons maintenant explorer le comportement de la sphere isotherme 
dans deux cas limites. 
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1.5 Cas limites 

Nous allons voir la limite A — ► qui correspond aux basses densites et 
aux hautes temperatures puis la limite A — > oo qui correspond aux hautes 
densites. 



1.5.1 Developpement a haute temperature 

Developpons la solution de l'equation de Lane-Emden ( 11.121 ) pres de 
A = en utilisant les conditions initiales ( 11.131) et ( 11.141) . On trouve 

^ ^—l x2 + w xi + jk- x6 + o ^- (L29) 

D'apres l'equation ( 11.211) . le developpement du contraste qui est le rapport 
de la pression au centre de la sphere et de la pression sur la paroi de la 
sphere est 

c = 1 + I A2+ iIo A4 + ^- 

Le developpement pour A petit correspond aux contrastes proches de 1 par 
valeurs superieures, c'est a dire a des configurations d'equilibre de la sphere 
isotherme presque homogenes. D'apres l'equation ( 11.171) . le developpement 
du parametre ij R est 

En se referant a la definition de 77 ( 11.171) . le developpement pour A petit 
correspond aux hautes temperatures. En utilisant les equations ( 11.181) et 
( 11.291 ). on trouve le developpement de / 

On en deduit la condition limite de l'equation differentielle ( 11.191) 

f(v R = 0) = 1 

et le developpement de / en t] R 

f(ri R ) = 1 - \v R ~ ^(r/T + 0[(^) 3 ] . (1.31) 
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D'apres les equations ( 11.231 ). ( 11.301) et ( 11.311) . le developpement de l'energie 
E est 



E 

NT 



= |-^A 2 + 0(A 4 ) = 5_3^ + 0[(r/T] 



D'apres l'equation ( 11.201 ). la pression du gaz en un rayon q inferieur au 
rayon Q de la paroi est telle que 



m P{q) 
Tp 



= 1 




+ 




1.32) 



Pour les hautes temperatures (A — > 0) : E ~ |iVT et P ~ l'energie 
cinetique l'emporte largement sur l'energie potentielle autogravitante et la 
pression suit la loi des gaz parfaits ; le systeme se comporte done comme un 
gaz parfait homogene. Les corrections par rapport a la loi des gaz parfaits 
dans l'equation ( 11.321 ) sont negatives puisque la force gravitationnelle est 
attractive. 

Nous allons maintenant explorer la limite A — ► oo de la sphere isotherme 
qui correspond a la limite des contrastes infinis. Cette limite approche une 
solution singuliere de l'equation de Lane-Emden isotherme ( 11.121) . 



1.5.2 Developpement autour de la sphere singuliere 

Presentons la limite A — > oo de la sphere isotherme [Sj. L'equation de 
Lane-Emden ( 11.121 ) a une solution singuliere 

2 

Xs(X) = In — . 

Cette solution diverge pour A = 0, elle ne verifie done pas les conditions aux 
limites ( 11.131 ) et ( 11.141) . On appelle cette solution de l'equation de Lane- 
Emden ( 11.121 ) solution singuliere et la sphere isotherme correspondante 
sphere singuliere. D'apres les equations ( 11.171 ) et ( 11.181 ). les valeurs des 
parametres rj R et / sont pour la sphere singuliere 

rj R = r]s = 2, f = f s = -. 

Etudions maintenant le comportement des solutions de l'equation ( 11.121) 
au voisinage de la solution singuliere. En posant 
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X = Xs + z 



la fonction z verifie l'equation 



fz 
dt 2 



dz 
dt 



2e 2 







Comme z <C 1, on a 



d 2 z dz 



dt 2 dt 
La solution de cette equation est 



-+2^=0 



z = A e 2 cos 



V7 



t + 5 



ou A et 5 sont des constantes d'integration. On a done 



X(A) = ln-^ + cos ( In A 



5 



La solution % approche de la solution singuliere xs pour A 
utilisant la relation ( 11.91 ). on trouve que la densite est 



oo. En 



Pm = Po tt; exp 



2 

A 2 " 



— r cos in A 



(5 



Pour A — > oo la densite se developpe en 



Pm Po 



A 2 



t4 / \/7 

1 + — r cos In A 

Al I 2 



5 



(A 



oo, 



D'apres cette equation et l'equation ( 11.211 ). on voit que la limite A — > oo 
correspond aux configurations d'equilibre de la sphere isotherme dont le 
contraste tend vers l'infini. On trouve que les parametres rj R ( 11.171) et / 
( 11.181 ) se developpent pour A — > oo en 
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7/ 



R 



= 2 



1 + 



A 



4Ai 



cos I —J- In A — S 



+ \[l sin ( In A — <) 



et 



/ = 



1 

3 



1 + 



,4 
4Al 

-V7 



3 cos 



V7 



In A - S 



sin 



'V7 



In A - 5 



On voit que la courbe f(r] R ) s'enroule en spirale autour du point singulier 

r)s = 2, fs = l 3 (figE|. 

Nous allons maintenant etudier la stabilite du gaz autogravitant a l'equi- 
libre thermodynamique. Notons tout de suite que dans la limite A — > oo 
etudiee dans ce paragraphe, le gaz autogravitant est instable. 



1.6 Stabilite 

Nous avons determine les configurations d'equilibre hydrostatique du 
gaz autogravitant a l'equilibre thermodynamique. Elles sont solutions de 
1' equation ( 11.111 ) et constituent l'approche du champ moyen de la meca- 
nique statistique des systemes autogravitants. Nous allons presenter dans 
cette section l'etude de la stabilite de ces configurations d'equilibre P, El, 
[IH[T2j. Lorsque ces configurations sont instables, le systeme cesse d'obeir 
a 1' equation ( 11.111) ; les particules s'effondrent sous l'effet de l'autogravite 
et collapsent en un point de densite infinie. Cette etude a mis en evidence 
que les regions de stabilite des configurations d'equilibre sont differentes 
suivant l'ensemble statistique dans lequel on se place. Nous nous limitons 
ici a l'ensemble microcanonique ou le systeme est isole thermiquement et 
a l'ensemble canonique ou le systeme est place dans un bain thermique. 
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1.6.1 Ensemble microcanonique et ensemble canonique 

L'ensemble microcanonique et l'ensemble canonique ne donnent pas le 
meme resultat en ce qui concerne la stabilite des configurations d'equilibre. 
Ceci s'explique par le comportement de grandeurs physiques comme les 
chaleurs specifiques et les compressibilites. 

Chaleurs specifiques des gaz autogravitants 

Les gaz autogravitants a l'equilibre thermodynamique ont des configu- 
rations d'equilibre de chaleur specifique negative (fig J1.2l dans le cas de la 
symetrie spherique). Cette propriete implique que le gaz a un comporte- 
ment inhabituel pour les systemes terrestres, il devient plus chaud quand 
il perd de l'energie. 

Dans l'ensemble canonique, ce comportement est source d'instabilite car 
le gaz echange de l'energie avec le thermostat. Supposons que la tempera- 
ture du gaz soit plus grande que la temperature du thermostat, un transfert 
d'energie s'opere du gaz vers le thermostat. Comme la chaleur specifique 
du gaz est negative, la temperature du gaz ne cesse de monter et l'ecart de 
temperature entre le gaz et le thermostat ne cesse de s'accroitre. Aucun etat 
d'equilibre ne peut etre atteint. Les configurations d'equilibre de chaleur 
specifique negative sont done instables dans l'ensemble canonique. Dans 
l'ensemble microcanonique, la situation est differente car le gaz est isole 
thermiquement . 

Nous sommes habitues aux systemes de particules en interaction a courte 
portee qui sont homogenes a l'equilibre thermodynamique. Pour ces syste- 
mes, l'ensemble canonique et l'ensemble microcanonique donnent les me- 
mes resultats dans la limite thermodynamique standard ou le nombre de 
particules N et le volume V verifient N — > oo, V — > oo avec y fini. 
L'energie de ces systemes est extensive (elle est proportionnelle au nombre 
de particules N et done au volume V dans la limite thermodynamique). 
On peut les diviser en sous-parties ; chaque sous-partie est en contact 
thermique avec le reste du systeme agissant comme un thermostat. Elle 
peut done etre consideree comme etant dans l'ensemble canonique, meme 
si le systeme est isole par rapport a l'exterieur. Les ensembles canon- 
ique et microcanonique sont done equivalents pour ces systemes et leur 
chaleur specifique doit etre positive. Les systemes autogravitants sont des 
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systemes de particules en interaction a longue portee qui ne sont pas ho- 
mogenes a l'equilibre thermodynamique. La limite thermodynamique stan- 
dard N —> oo, V — > oo avec y fini n'est pas pertinente. La limite thermody- 
namique pertinente des systemes autogravitants est N — ► oo, V — > oo avec 
-^r qui est fini. Leur energie est extensive, elle est proportionnelle au nom- 
bre de particules N mais n'est plus proportionnelle au volume V dans la 
limite thermodynamique autogravitante. L'ensemble canonique et l'ensem- 
ble microcanonique ne donnent pas les memes resultats en ce qui concerne 
la stabilite des configurations d'equilibre des gaz autogravitants ; ils ont 
des configurations d'equilibre stable avec une chaleur specifique negative 
dans l'ensemble microcanonique fTTl fT^j. 

En revanche, pour tous les systemes autogravitants ou non la chaleur 
specifique d'un systeme est toujours positive dans l'ensemble canonique. 
Soit un systeme place dans un bain thermique a la temperature T et ayant 
des niveaux d'energie E{. Son energie moyenne est 

<E>= £iflexp(-/?£Q 
£ ? exp(-/^) 

ou (3 = ^ est le facteur de Boltzmann. Sa chaleur specifique est egale au 
carre des fluctuations de l'energie. En effet, 

d < E > 2 d< E> 2 2 
C = ^^ = - ^^ = /? <(E~<E>f> . 

La chaleur specifique d'un systeme est toujours positive dans l'ensemble 
canonique. Une telle contrainte sur le signe de la chaleur specifique n'existe 
pas dans l'ensemble microcanonique. 

Nous allons maintenant voir les proprietes des compressibilites isotherme 
et adiabatique. 



Compressibilites des gaz autogravitants 

La compressibilite d'un fluide mesure le rapport entre sa variation de 
volume V et sa variation de pression P. Son expression est 

l_dV 
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Lorsque k > 0, on a |p < 0, c'est a dire que le fluide diminue de volume 
quand on le comprime plus, un tel comportement est normal pour un fluide. 
Par contre lorsque k < 0, on a |p > 0, c'est a dire que le fluide diminue 
de volume quand on le comprime moins, un tel comportement est anormal 
pour un fluide et il conduit a une instabilite analogue a celle de Jeans j2Bj 
pour un fluide autogravitant homogene (voir annexe B). Ainsi l'etude du 
signe de k permet de savoir si les configurations d'equilibre des systemes 
autogravitants sont des configurations d'equilibre stables ou instables. Le 
systeme est stable lorsque k > et instable lorsque k < 0. 

Dans l'ensemble canonique ou la temperature est constante, c'est la 
compressibilte isotherme qui doit etre positive et dans l'ensemble micro- 
canonique ou l'entropie est constante, c'est la compressibilte adiabatique 
qui doit etre positive. 

Nous allons determiner les configurations d'equilibre stable de la sphere 
isotherme dans ces deux ensembles. 

1.6.2 Stabilite de la sphere isotherme 

Chaque configuration d'equilibre de la sphere isotherme est caracterisee 
par une valeur du rayon reduit A ( 11.101) . et a chaque valeur de A correspond 
une valeur du contraste C ( 11.211) qui est le rapport de la pression au centre 
de la sphere et de la pression a la peripheric de la sphere. La limite A — > 
correspond au gaz parfait homogene pour laquelle C — > 1, il s'agit de la 
limite des basses densites. La limite A — > oo correspond a la sphere sin- 
guliere pour laquelle C — > oo, il s'agit de la limite des hautes densites. Nous 
allons determiner la zone de stabilite de la sphere isotherme, c'est a dire les 
valeurs du rayon reduit A et celles du contraste pour lesquelles les configu- 
rations d'equilibre sont stables et determiner les points qui correspondent 
a ces valeurs dans le diagramme de phase (fig Jl.lj) . 

Determinons d'abord la zone de stabilite dans l'ensemble canonique. 

Stabilite de la sphere isotherme dans l'ensemble canonique 

La zone de stabilite dans l'ensemble canonique correspond aux valeurs 
du rayon reduit A comprises entre et 6.5..., c'est a dire aux valeurs du 
contraste C comprises entre 1 et 14.1.... Ces configurations d'equilibre sta- 
ble correspondent a la partie superieure du diagramme de phase (fig Jl.ll) 
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du point (rj R = 0,/ = 1) au point (rj R = 2.43...,/ = 0.40...). Ces configu- 
rations sont stables dans 1' ensemble canonique car les chaleurs specifiques 
et la compressibilite isotherme sont positives (fi g)1.2j) . En revanche, pour 
A = \ can = 6.5... (C = 14.1..., rf = j] can = 2.43...,/ = 0.40...), la chaleur 
specifique a pression constante et la compressibilite isotherme divergent 
en devenant negatives. Cette configuration est le point d'instabilite de la 
sphere isotherme dans l'ensemble canonique. D'apres les equations ( 11.171) . 
( 11.181) et ( 11.271) . les equations definissant r] can et X can sont 

6 f{j] mn ) = r] can . (1.33) 

et 

Mean) _ (* ( A ^)) _ n 

2 

Determinons maintenant la zone de stabilite dans l'ensemble micro- 
canonique. 

Stabilite de la sphere isotherme dans l'ensemble microcanonique 

La zone de stabilite dans l'ensemble microcanonique correspond aux 
valeurs du rayon reduit A comprises entre et 25.8..., e'est a dire aux valeurs 
du contraste C comprises entre 1 et 389. Ces configurations d'equilibre sta- 
ble correspondent a toute la partie superieure du diagramme de phase 
(figO) du point {r] R = 0,/ = 1) au point (r) R = 2.51.. .,/ = |) et 
a la inferieure du diagramme du point (77^ = 2.51..., / = |) au point 
(r] R = 2.14..., / = 0.26...). Ces configurations sont stables dans l'ensem- 
ble microcanonique car la compressibilite adiabatique est positive (fig jl.2l ). 
Pour A = X mic = 25.8 (C = 389, rj R = rj mic = 2.14...,/ = 0.26...), la com- 
pressibilite adiabatique diverge en devenant negative. Cette configuration 
est le point d'instabilite de la sphere isotherme dans l'ensemble micro- 
canonique. On peut remarquer aussi que la chaleur specifique a pression 
constante s'annule en redevenant positive. D'apres les equations ( 11.171 ). 
( 11.181) et ( 11.281) . les equations definissant 7/ m j c et X m i c sont 

48/(r/ TOJC ) 2 + (Sf] mic - 38)/(r/ TOiC ) + j] mm = . (1.34) 

et 
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48A^ C e 2 *^ + 3 [8A m , cX '(A m , c ) + 38] \ mic e x{Xmic) xX^mic 

9A TO j C [^ (A m ^ c )] 



Les configurations d'equilibre correspondant a la la partie superieure du 
diagramme de phase (fig JI.ll) du point (rj R = 2.43...,/ = 0.40...) au point 
(r) = 2.51..., / = g) et a la partie inferieure du diagramme du point (rj R = 
2.51..., / = |) au point (r] R = 2.14..., / = 0.26...) ont une chaleur specifique 
a pression constante negative. Elles sont stables dans dans l'ensemble mi- 
crocanonique et instables dans l'ensemble canonique. La zone de stabilite 
dans l'ensemble microcanonique est plus etendue que dans l'ensemble canon- 
ique car le systeme y subit plus de contraintes. Ces resultats sur la stabilite 
de la sphere isotherme, deduits du comportement des chaleurs specifiques 
et des compressibilites sont confirmes par les calculs Monte Carlo faits dans 
l'ensemble canonique et dans l'ensemble microcanonique |6j. 



1.6.3 Instabilites gravitationnelles en astrophysique 

On voit done que l'ensemble statistique joue un role important dans la 
physique des gaz autogravitants. On doit done determiner dans quel en- 
semble statistique se trouvent les objets astrophysiques que Ton etudie. 
Les nuages interstellaires et les distributions de galaxies sont baignes par 
le rayonnement de fond micro-onde qui joue le role de thermostat. On doit 
done utiliser l'ensemble canonique pour les etudier. Les etoiles sont des 
systemes isoles, on doit done utiliser l'ensemble microcanonique pour les 
etudier. Les instabilites gravitationnelles permettent d'expliquer la forma- 
tion des etoiles dans les nuages interstellaires [2HI EDI • Les instabilites grav- 
itationnelles jouent aussi un role important dans la formation de structures 
hierarchiques dans les distributions de galaxies, les amas et les superamas 
de galaxies. 
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Chapitre 2 

Mecanique statistique 



La mecanique statistique des systemes autogravitants est l'objet de ce 
chapitre. Nous nous placerons dans le formalisme de Gibbs pH El EH E| 
plutot que dans le formalisme de Boltzmann (TTJ H2] • Elle a ete etudiee dans 
l'ensemble microcanonique, dans l'ensemble canonique et dans l'ensemble 
grand-canonique [EH El IE|- La mecanique statistique dans l'approche du 
champ moyen decrit exactement la phase gazeuse dans la limite thermo- 
dynamique autogravitante ou le nombre de particules N et le volume V 
tendent vers l'infini et ou est fini. Elle montre que le gaz autogravitant 

obeit a l'equation d'equilibre hydrostatique ( 11.111) et obeit localement a 
1' equation d'etat des gaz parfaits. Nous allons presenter la mecanique statis- 
tique dans l'ensemble microcanonique et dans l'ensemble canonique. Dans 
la limite thermodynamique autogravitante, la fonction de partition est ap- 
prochee par une integrale fonctionnelle sur la densite. Le poids statistique 
de chaque densite est l'exponentielle d'une "action effective" proportion- 
nelle a N. Dans la limite N — > oo, on applique 1' approximation de point col 
qui constitue l'approche du champ moyen. Les points col sont identiques 
dans l'ensemble microcanonique et dans l'ensemble canonique. lis corre- 
spondent aux configurations d'equilibre hydrostatique des systemes auto- 
gravitants isothermes dont la densite de masse obeit a l'equation ( 11.111 ). 
Si le poids statistique de la densite dans l'integrale fonctionnelle diminue 
pour de petites fluctuations autour du point col, alors le point col domine 
l'integrale et le champ moyen est valide dans la limite iV — ► oo. Par contre, 
si le poids statistique de la densite dans l'integrale fonctionnelle augmente 
pour de petites fluctuations autour du point col alors le point col ne domine 
pas l'integrale et le champ moyen n'est pas valide dans la limite iV — > oo. 
Les calculs Monte Carlo jHj confirment les hypotheses du premier chapitre 
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sur la stabilite des configurations d'equilibre. Dans l'ensemble canonique, 
le champ moyen cesse d'etre valide lorsque la compressibilite isotherme 
diverge et devient negative. Dans l'ensemble microcanonique, le champ 
moyen cesse d'etre valide lorsque la compressibilite adiabatique diverge et 
devient negative. Dans ces conditions, il faut etudier le systeme par des 
calculs Monte Carlo. Pour ces points d'instabilite, les calculs Monte Carlo 
montrent qu'une transition de phase se produit de la phase gazeuse vers la 
phase collapsee. 



2.1 Ensemble microcanonique 

Etudions la mecanique statistique dans l'ensemble microcanonique d'un 
gaz thermiquement isole dans un volume V et d'energie E, compose de N 
particules de masse m interagissant entre elles par la gravite. Calculons 
d'abord le nombre de microetats. 



2.1.1 Nombre de microetats 

Le hamiltonien du gaz dont les particules ont comme positions q\, 
et comme impulsions pi, est 



N -2 i 

(2-1) 

ou le premier terme est l'energie cinetique et le deuxieme terme Ep est 
l'energie potentielle. Comme les forces non gravitationnelles dominent a 
courte distance dans les systemes physiques qui nous interessent (nuages 
interstellaires, distributions de galaxies), on a introduit le cut-off suivant 



Wi-QjU = Wi-Qj\ si \qi-qj\>A 



et 



Wi-qj\A = -A si \qi-qj\<A. (2.2) 
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Le nombre de microetats en fonction de l'energie E est 



N 

n 

1=1 



(2tt) 3 



6 



E 



N 

E 

i=i 



Pi 

2 m 



- E P (q h ...,q N ) 



(2.3) 

Par convention, la constante de Planck h est choisie egale a 1. Notez que 
le cut-off a courte distance A permet d'eviter la divergence de Q(E) qui 
ne serait pas defini mathematiquement. En effectuant le changement de 
variable polaire sur les impulsions p 2 = Yli=i <T™ e ^ en integrant sur les 
angles, on trouve 



On introduit les positions sans dimension 



l'energie potentielle sans dimension 



u(r h ...,r N ) = --^—E P (q L ,...,q N ) = 1 £ 1 (2.4) 



1XV v „ Pi ' 7 1 a 
1<Kj<N ' J ' 



avec a = -4- <C 1 et le parametre 

^3 



Gm 2 N 2 



(2.5) 



On trouve que le nombre de microetats ft(E) est egal au produit du nom- 
bre de microetats Qqp du gaz parfait (GP) d'energie E et de volume V 
contenant N particules de masse m par une integrale sur les positions des 
particules Qi nt qui contient l'information sur l'interaction gravitationnelle, 
ainsi : 
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Q(E) = Qgp fynt > 

~ iVT ~Y~ \YJ V E ' 

N - — 1 

n* = £ ¥-./ ndS J f+ M) ! ^ (e + l). (2 . 6) 



/=1 



La presence de la fonction de Heavyside 6(e + 4?) impose que l'energie 
potentielle doit etre inferieure a l'energie totale E. Connaissant le nombre 
de microetats (eq. (12.61) ). nous en deduisons l'entropie S = \nft(E) et toutes 
les grandeurs physiques. 



2.1.2 Grandeurs physiques 

L'entropie du systeme est 



S = lntt(E) = lntt GP + Inn 



int 



(2.7) 



Elle s'exprime comme la somme de deux termes. Le premier terme Sgp = 
IuQgp est l'entropie d'un gaz parfait d'energie E et de volume V compose 
de N particules. Dans la limite N — > oo, on obtient la formule de Sackur- 
Tetrode EH 



S GP = N 



, V 3, fmE\ 5 

In h - In + - 

N 2 \37riV/ 2 



Le deuxieme terme de l'equation (12.71) contient l'information sur l'interac- 
tion gravitationnelle entre les particules. 

A partir de l'entropie, on calcule la temperature T et la pression de paroi 
P par les relations thermodynamiques standard ^ = (f§)y )iV et P = 
T (§^) EN - En utilisant les equations (12.51) . (12.61) et (12.71) . on trouve 



1 _ 3 TV ed SN 

t ~ Ye + Ed~e { mt) ~ Ye 



3N 



1 



N 



e+ 



E 



ou 
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/ \ UL _ 9 

/ 1 \ = koi^^e n £id 3 n (e + ft) ' - + 

\ e + # / /„ oi „ m e „„«e n £id 3 n (e + ft) ¥ - 1 «(€ + ft) ' 

Grace aux equations ( 12.51) . ( 12.61) et ( 12.71 ). on trouve l'equation d'etat 



= 1 + (In Q in *) = - + + - 

NT 3 N de K mt) 2 3 N 2 

(2-9) 

En combinant les equations ( 12.81) et ( 12.91 ). on retrouve le theoreme du viriel 

(H22D 

3PV = -NT + £7. 
2 

Nous allons maintenant nous placer dans l'ensemble canonique. 

2.2 Ensemble canonique 

Etudions la mecanique statistique dans l'ensemble canonique d'un sys- 
teme de volume V compose de N particules de masse m interagissant entre 
elles par la gravite. II est place dans un bain thermique a la temperature 
T et une pression P s'applique sur la paroi qui l'enferme. 




2.2.1 Fonction de partition 

La fonction de partition est 

* - ^ / - ("0) 

j i=i v ' 

ou H est le hamiltonien du gaz defini dans la section precedente par 
l'equation ( 12.11 ). Remarquons que le cut-off a courte distance dans l'energie 
potentielle permet de definir mathematiquement les integrales dans Z . En 
calculant les integrales gaussiennes sur les impulsions et en introduisant les 
positions sans dimension des particules f\ = -pr, on trouve que la fonction 

de partition Z est egale a la fonction de partition Zqp du gaz parfait (GP) 
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de temperature T et de volume V contenant N particules de masse m fois 
une integrate sur les positions des particules Zi nt qui contient l'information 
sur l'interaction gravitationnelle, c'est a dire precisemment : 



Z — Zgp Z in t , 



N 



z, 



int 



'volume unite 



J~J d 3 r/ e V u(r r N ) 



Le parametre 77 vaut 



77 = 



Gm 2 N 
T V% 



[2.12) 



et u est l'energie potentielle sans dimension definie par l'equation ( 12.41) . 
Connaissant la fonction de partition Z (eq.f l2.llT )). nous allons maintenant 
en deduire l'energie libre F = —TlnZ et toutes les grandeurs physiques. 

2.2.2 Grandeurs physiques 

L'energie libre F = —TlnZ est d'apres l'equation ( 12.111 ) 



F = F GP - T $ n (ti) 



[2.13) 



ou Fgp = — T In Zqp est l'energie libre du gaz parfait de temperature T et 
de volume V contenant TV particules de masse m. Dans la limite N — > 00, 
on obtient 



F GP = -NT In 



eV ( mT 



N V 2 7T 



D'apres les equations ( 12.111 ). ( 12.121 ) et ( 12.131) . la pression P 
qui s'exerce sur le systeme est 



NT 77 T / 



(dF\ 
\dV/T,N 



;2.i4) 
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D'apres les equations ( 12.111) et ( 12.121) . l'energie moyenne < E >= 

3 



< E > 



NT — T r] <5> N ( V ) 



;2.i5) 



En combinant les equations ( 12.141 ) et ( 12.151) . on retrouve le theoreme du 
viriel (TT32T) 



3PV=-NT+ < E> 
2 



Introduisons la grandeur / 



PV 



1 



(2.16) 



NT 3N 

II s'agit de l'equation d'etat du systeme. Dans la limite i] — > (gaz parfait), 
on a f(r] = 0) = 1 et <&jv(?7 = 0) = 0. En integrant la relation ( 12.161) . on a 



<$> N ( V ) = 3 N P dx - — ^ 
Jo x 



[2.17) 



On peut exprimer toutes les grandeurs thermodynamiques en fonction de 
la grandeur f(rj) grace a cette equation. D'apres l'equation ( 12.131 ). l'energie 
libre est 



F = Fqp 



3NT 



'' dx l^M 



Jo x 

L'energie moyenne est en utilisant l'equation ( 12.151 ) 



;2.i8) 



< E > = 3 NT 



(2.19) 



On en deduit la valeur de l'entropie S = 



E-F 
T 



S = S GP + 3N 



m - 1 + [ v dx — M 

JO x 



(2.20) 



Calculons la chaleur specifique a volume constant et la chaleur specifique a 
pression constante. D'apres l'equation ( 12.201) . la chaleur specifique a volume 
constant c v = (^) VN vaut 



Cii — 3 



fin) - vf'W - i 



12-21) 
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D'apres les equations ( 12.121 ). ( 12.141 ) et ( 12.211 ). la chaleur specifique a pression 
constante [27J 



J- \dTJV,N 



\dV/T,N 



devient 



p ~ * m + Win) ' ( ' 

On peut aussi calculer la compressibilite isotherme et la compressibilite 
adiabatique. En utilisant les equations ( 12.121 ) et ( 12.141) . la compressibilite 
isotherme kt = — y {jfp) T N s'exprime suivant 

NT 1 

~ KT = m+bm ■ (2 ' 23) 

La compressibilite adiabatique 



1 fdV\ Cy 

S,N C P 



= 77 I 7777 J = — K T , 

V \ OF / c at C r 



d'apres les equations (E21) et (E2BD vaut 

^ Ks = *JM - - I] ( 2 24) 

Nous allons maintenant exposer l'approximation de champ moyen qui 
qui est exacte dans la limite N — > oo et qui, nous allons le voir, conduit a 
l'hydrostatique presentee dans le premier chapitre. 



2.3 Champ moyen 

Presentons le champ moyen dans l'ensemble canonique puis dans l'ensem- 
ble microcanonique. 
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2.3.1 L'ensemble canonique 

On va se placer dans la limite N — > oo. Nous allons montrer que dans 
cette limite, la fonction de partition ( 12.111) devient une integrale fonction- 
nelle sur la densite dont le poids statistique de chaque densite est l'ex- 
ponentielle d'une "action effective" proportionnelle a N [19J, en utilisant 
l'approche exposee dans ref.[32j. Pour cela, on va diviser le volume unite 
d'integration de la fonction de partition ( 12.111) en M cellules (1 « M « iV) 
de volume jy suffisamment grandes pour contenir un grand nombre de par- 
ticules et suffisamment petites pour que le potentiel gravitationnel puisse 
etre considere comme uniforme dans chaque cellule, reintegration sur les 
positions f\, r)v devient une somme discrete sur le nombre de particules 
Til, ■ ■ • , tim par cellule. On a 



ni,n 2 ,...,n M 



ni!n 2 !...nM! v ^ ' \M J 



x exp — n a u ab n b 



N 

a.b 



OU 



Uab 



i - n 



a 



f a et fb etant la position respective du centre de la cellule a et du centre de 
la cellule b. Le volume d'integration elementaire Yii=i d 3 n devient ( -p) 
et — r- — r est le nombre de combinaisons qui ne sont pas equivalentes 
pour placer les Tii, Ti2,...,n m particules dans chaque cellule. En utilisant la 
formule de Stirling n\ ~ y/2 tt n n n e~ n pour n —> oo, on a 



iV>l 



ni,n 2 ,...,n M 



x exp 



y^n a In 



'n a M' 
N 



+ — £^n a u ah n h 



2N 



Introduisons la densite de particules p(r) qui vaut 7Jl jf- sur la cellule a. On 
obtient 
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N(l - / d 3 rp(f)) 



d& 

2^ eXP 



iNb( J dVp(r) - r 



L'integrale e® N ^ se transfome en une integrate fonctionnelle sur la densite 
de particules p(r) 



iV>l 



/A 1^ 
Dp(.)— exp[-N s c (p(.),b)] (2.25) 



avec 1' "action effective" 



s c (p(.),b) = J d¥p(r) lnp(r) - 

— y d 3 f /o(r^ 



77 / d r d r 



7T p(0 



J2.26) 



L'integrale fonctionnelle ( 12.251) est dominee pour N — ► 00 par le point 
col de 1' "action effective" ( 12.261) qui verifle les relations suivantes 



ds f 



(pcoh hoi) = 



5s r 



{pcoh hoi) = 



La premiere relation impose la normalisation de la densite 



J d 3 r Pcol {r) = 1 



[2.27) 



La seconde relation impose que la densite soit solution de l'equation de 
point col 



\np cd (r) - 77 



777 Pcoi{r ) = a co i 



[2.28) 



a C oi = i b co i — 1 est un multiplicateur de Lagrange associe a la condition de 
normalisation de la densite ( 12.271 ). En appliquant le Laplacien a l'equation 
de point col et en introduisant la fonction $(f) = In p co i(f), on trouve 



V^>(r) + 4 77 77 e $(f) = 



2.29) 
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II s'agit de l'equation de Liouville avec un signe oppose par rapport a 
l'equation en theorie des champs sans gravite. Le signe ici correspond a 
l'attraction de la force gravitationnelle et la theorie avec l'instabilite est 
le secteur physique. La densite sans dimension p co i{r) = est liee a la 
densite de masse p m (q) introduite dans le chapitre 1 par les relations 

p m (q) = — Pcolif) = — e m > q = V*r. (2.30) 

En utilisant les equations ( 12.121) et ( 12.301) . on trouve que l'equation de point 
col est identique a l'equation d'equilibre hydrostatique ( 11.111) 

t42 i 4 7T G m 

Vf\np m + p m {q) = 0. 

Les solutions de point col sont done equivalentes aux configurations d'equi- 
libre hydrostatique du systeme autogravitant isotherme obeissant a l'equa- 
tion des gaz parfaits. Rappelons qu'en hydrostatique, on devait supposer 
l'equation d'etat. Dans la limite thermodynamique diluee ou N — > oo, 
V — > oo et -p- est fini, la mecanique statistique demontre que le gaz obeit a 

l'equation d'etat des gaz parfaits (CD) ; ce resultat a ete determine grace aux 
informations microscopiques de la mecanique statistique donnees par les 
forces de Newton s'exergant entre les particules. D'apres l'equation ( 12.301) . 
la pression au point q = V$r est 

P(q) = ■ ( 2 - 31 ) 

Dans le cas de la symetrie spherique, la grandeur / = jtL (eq. ( 12.161) ) 
coincide avec la grandeur introduite dans l'equation ( 11.181 ) introduite dans 
le cadre de l'hydrostatique et obeissant a l'equation differentielle ( 11.191) . 
On a 

NT J 

L'equation ( 11.191) s'integre de cette maniere 

3 / - (1 - f(x)) = 3 [f(ri R ) - 1] +7? - ln/(^) . 
Jo x 
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En utilisant les equations ( 12.181) et ( 12.201) . on en deduit que l'energie libre 
et l'entropie verifient les relations suivantes 



= 3 [1 - f(r) R )] - V R + In f( V R ) , (2.32) 

= 6 [f( V R ) - 1] + n R - In /(^) . (2.33) 

En utilisant les equations (l2~m (l2~m (E22), (E2D et ( 12~M ). on retrouve 
les expressions de l'energie, de la chaleur specifique a volume constant, de 
la chaleur specifique a pression constante, de la compressibilite isotherme 
et de la compressibilite adiabatique des configurations d'equilibre hydro- 
statique. 

Nous allons maintenant presenter l'approche du champ moyen dans 
l'ensemble microcanonique. 

2.3.2 L'ensemble microcanonique 

Exprimons l'integrale sur les positions Qi n t dans le nombre de microetats 
defini par l'equation ( 12.61) en terme de l'integrale e® N ^ definie par l'equa- 
tion ( 12.111) . Pour cela, on utilise la transformation de Fourier suivante j23| 



= ^ r^T^I- (2-34) 
D'apres les equations (ESI), (EUD) et ( l2~m on a 

i 

ou 6 = (Fm^N 2 (eq- (l2.5l) ). On introduit la variable d'integration r] = ^ 



7 



ou le contour d'integration 7 est un contour parallele a l'axe des imaginaires 
purs. En utilisant la formule de Stirling pour la fonction T, on trouve que 
pour N ^> 1 

Q int = I ^. e ^ + *»W-W). (2.35) 

7 7 27T2 
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En utilisant la representation de l'integrale sur les positions e $n ™ par 
l'integrale fonctionnelle ( 12.251) . on obtient 

Q mt = f Dp(.) — pr exp[-N s mic (p(.),b,ri)} (2.36) 
J Ztt Zm 

avec 1' "action effective" 



i _ >x 1 I ^ T] f d 3 f d 3 f' 



Smic(p{-),b,r]) = J dfp(r) lnp(r) - - j p(r) p(f ) 

+ ib [1 - J d 3 f p(r)] + ^ In?/ - rj e . (2.37) 

On a transformed le nombre de microetats en une integrale fonctionnelle 
sur la densite. L'integration sur b contraint la normalisation de la densite 
comme dans la fonction de partition ( 12.251 ). L'integration sur j] contraint 
l'energie comme cela doit etre le cas dans l'ensemble microcanonique. 

Comme pour la fonction de partition ( 12.251 ). on va appliquer au nombre 
de microetats ( 12.361) l'approximation de point col. L'integrale fonctionnelle 
( 12.361) est dominee pour TV — > oo par les extrema de 1' "action effective" 
( 12.371 ) qui verifient les relations suivantes 

filmic q filmic q d^mic q 

db 5p(-) drj 

Comme dans l'ensemble canonique, la premiere relation impose la normal- 
isation de la densite ( 12.271 ) et la deuxieme relation impose que la densite 
soit solution de l'equation de point col ( 12.281) . La troisieme relation impose 
la contrainte suivante sur r] 

3 If d 3 f d 3 f ' , . , ./. 

7T Pif) P(r ) • (2.38) 




2t] 2 J \r — r 
En utilisant les equations (12.51 ). ( 12.121) et ( 12.301 ). cette relation devient 

3^ G fd 3 qd 3 q 
E = -NT + — / — — — p m (q) p m {q ) • 
2 2 J \q- q \ 

Ainsi la relation ( 12.381 ) contraint l'energie dans l'ensemble microcanon- 
ique. Les points col sont les memes dans l'ensemble microcanonique et 
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dans l'ensemble canonique et verifient la condition d'equilibre hydrosta- 
tique ( II. llj ). L' approximation de champ moyen dans 1'ensemble canonique 
et l'approximation de champ moyen dans l'ensemble microcanonique cor- 
respondent done aux configurations d'equilibre de l'hydrostatique quand 
N -> oo . 

2.4 Calculs Monte Carlo 

L'algorithme de Metropolis JMJ a ete applique aux systemes autograv- 
itants isothermes dans un cube de volume V dans l'ensemble canonique 
et dans l'ensemble microcanonique, pour un nombre de particules allant 
jusqu'a 2000 Un cut-off a courte distance a (a ~ 10~ 3 — 10~ 6 ) a ete 
introduit dans l'interaction gravitationnelle. 

2.4.1 Algorithme de Metropolis 

L'algorithme de Metropolis permet de simuler la thermalisation d'un 
systeme et de calculer les moyennes thermodynamiques de ces grandeurs 
physiques. Placons-nous dans l'ensemble canonique, le systeme etant en 
contact avec un thermostat a la temperature T. Pour plus d'explications, on 
peut se referer au livre de Binder et Heerman . La fonction de partition 
Z du systeme (eq. ( 12. lip ) est le produit de la fonction de partition du gaz 
parfait Zqp et de l'integrale sur les coordonnees des particules e® N ^ 

N 

Z = Z GP e* N W , e *»W = / JJ d 3^ e 77<n,...,f w ) _ 

J volume unite 

A l'equilibre thermodynamique, chaque configuration du systeme x = 
(fi, ...,r/v) a une probabilite d'avoir lieu, qui est 

p rju{x) p—T- 
P ( x ) = - = 

r ^y x ) e <& N ( n ) e $ N { n ) 

ou Ep = —Tt] u(x) est l'energie potentielle de la configuration x du 
systeme. Expliquons maintenant comment on simule la thermalisation du 
systeme a partir de multiples transformations d'une configuration initiale 
prise au hasard. Une transformation de la configuration x = (fi, r/v) con- 
siste a changer au hasard la position de l'une des particules. Soit W(x — > x) 
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la probability de transition d'une transformation d'une configuration x vers 
une nouvelle configuration x . Si 



P eq {x) W(x -> x) = P eq (x) W(x -> x) 

alors la probabilite de la configuration x tend vers la probability d'equilibre 
P eq (x) apres un grand nombre de transformations. On doit done avoir 



w ( x Zl e vSu = 

W(x'^x) 



ou 6Ep = —Trj 6u = —Tr] (u(x') — u{x)) est la variation d'energie poten- 
tielle entre les deux configurations x et x . On peut choisir par exemple 

W(x^>x) = e 11 ^ si 6u<0 {8E > 0) 



W(x -+x) = l si 5u>0 (5E < 0). 

Presentons maintenant l'algorithme de Metropolis 

1. On selectionne au hasard un changement de configuration x — > x . 

2. On calcule la variation d'energie potentielle 5u correspondant a ce change- 
ment de configuration . 

3. On calcule la probabilite de transition W correspondant a ce changement 
de configuration . 

4. On tire au sort un nombre z entre et 1. 

5. Si z < W alors on effectue le changement de configuration x — ► x , sinon, 
on conserve la configuration x. 

Ce processus constitue un "tour" Monte Carlo. Au bout d'un grand nombre 
de "tours", le systeme est a l'equilibre thermodynamique. A partir de la, 
on effectue M tours supplementaires et on calcule la moyenne statistique 
d'une grandeur A grace a l'equation suivante 

1 M 

<A> = mH A ^ ' 

i=i 

chaque configuration xi intervenant avec la probabilite P eq (xi). 
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2.4.2 Resultats dans l'ensemble canonique 

Presentons les resultats des calculs Monte Carlo des systemes autograv- 
itants thermalises dans 1'ensemble canonique. 
II y a deux phases distinctes : 

-une phase gazeuse pour r] R < 2.43... 

-une phase colllapsee pour r] R > 2.43.... 
Dans la phase gazeuse, les calculs Monte Carlo sont insensibles au cut-off 
a courte distance de la gravite et reproduisent remarquablement bien les 
resultats du champ moyen a partir d'un nombre de particules assez bas 
(N > 200). Pour f] R ~ 2.43, il y a une brutale transition de phase. La pres- 
sion devient grande et negative. Les particules sont aspirees vers le centre 
et tout le systeme collapse en un corps tres dense. Le point de collapse 
77 ~ 2.43 obtenu par les calculs Monte Carlo correspond au point d'in- 
stabilite r] can = 2.43... (eq. ( 11.331) ) de 1' ensemble canonique prevu dans le 
chapitre 1 ou la chaleur specifique a pression constante c p et la compress- 
ibilte isotherme kt divergent. 

2.4.3 Resultats dans l'ensemble microcanonique 

Presentons maintenant les resultats des calculs Monte Carlo des syste- 
mes autogravitants thermalises dans l'ensemble microcanonique. Comme 
dans l'ensemble canonique, les calculs Monte Carlo reproduisent remar- 
quablement bien les resultats du champ moyen a partir d'un nombre de 
particules N > 200. Le point de collapse obtenu par les calculs Monte 
Carlo est le point 77 ~ 2.17 dans la partie inferieure du diagramme de 
phase f(r) R ) = (fig Jl.il) . Ce point est tres proche du point d'insta- 
bilite de l'ensemble microcanonique j] m i c = 2.14... (eq. ( 11.341) ) prevu dans 
le chapitre 1 ou la chaleur specifique a pression constante c p s'annule et la 
compressibilite adiabatique k$ diverge. II se produit a ce point la transi- 
tion de la phase gazeuse vers la phase collapsee et les calculs Monte Carlo 
montrent que la phase collapsee prend la forme d'un coeur-halo ou les par- 
ticules sont tres condensees au centre et ou un halo de particules demeure 
a la peripheric du systeme. Dans l'ensemble canonique, la phase collapsee 
a la forme d'un corps tres dense ou sont concentrees toutes les particules. 
La phase collapsee prend done des formes differentes dans l'ensemble mi- 
crocanonique et dans l'ensemble canonique. 



48 



Les calculs Monte Carlo reproduisent remarquablement bien les resul- 
tats du champ moyen et permettent d'etudier la phase collapsee. Tous les 
resultats dans la phase gazeuse sont insensibles au cut-off a courte distance. 
Le calcul de la fonction de partition est moins sensible a la singularity a 
courte distance en de l'interaction gravitationnelle que la resolution des 
equations du mouvement de Newton pour N particules. La determination 
du mouvement classique de N particules interagissant par la gravitation 
ou la resolution des equations de Boltzman comprenant l'interaction grav- 
itationnelle a N particules demandent des algorithmes sophistiques pour 
eviter de trop longs calculs numeriques |36j ; elles donnent des informa- 
tions sur le comportement dynamique en dehors de l'equilibre thermique. 
Cependant les calculs Monte Carlo sont largement sumsants pour decrire 
les systemes autogravitants et en donnent une comprehension approfondie. 
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Chapitre 3 



Systemes autogravitants comportant 
plusieurs sortes de particules 



Nous avons j usque la considere les systemes autogravitants composes 
par des particules ayant toutes la meme masse. Dans l'univers, les systemes 
que Ton peut considerer comme autogravitants sont souvent composes de 
particules ayant des masses differentes. C'est le cas, d'une part des nuages 
interstellaires composes de plusieurs types d'atomes et de molecules, c'est le 
cas, d'autre part des distributions de galaxies, les galaxies ayant des masses 
differentes. II est done important d'etudier les systemes autogravitants a 
l'equilibre thermodynamique composes de plusieurs sortes de particules, 
chaque sorte de particules ayant une masse differente. Nous allons etudier 
d'abord l'hydrostatique d'un melange de gaz autogravitants, chaque gaz 
autogravitant etant compose par des particules ayant la meme masse et 
obeissant localement a l'equation d'etat des gaz parfaits. Ensuite, nous al- 
lons etudier la mecanique statistique des systemes autogravitants composes 
de plusieurs sortes de particules et montrer que l'approximation de champ 
moyen conduit a l'hydrostatique d'un melange de gaz autogravitants et 
permet de deriver les equations d'etat du systeme [15j. 



3.1 Hydrostatique 

Considerons un systeme autogravitant isotherme constitue par un me- 
lange de gaz obeissant localement a l'equation d'etat des gaz parfaits. II est 
compose de gaz de Ni particules de masse rrii (1 < i < n). Etant donne la 
densite de masse p mi (q) des particules de masse m z au point q, la pression 
partielle des particules de masse m 2 au point q est (J3D 
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T 

Pi(q) = — pmM ■ 

Presentons maintenant la condition d'equilibre hydrostatique qui s'appli- 
que a ce melange de gaz autogravitants isothermes. 

3.1.1 Equilibre thermodynamique 

Nous allons exprimer l'equilibre de chacun des gaz de particules de masse 
Tfii (1 < % < n). Pour que les particules de masse m 2 soient en equilibre 
hydrostatique, il faut que les forces de pression s'appliquant sur elles com- 
pensent les forces gravitationnelles s'appliquant sur elles. Cette condition 
(partielle) d'equilibre hydrostatique s'ecrit 

"V^ + PmMM = (3.1) 

ou g(q) est le champ gravitationnel au point q engendre par l'ensemble du 
melange. En introduisant la pression totale 



p(q) = j2m (3.2; 



i=l 



et la densite totale 



Pm(q) = J2pmM ( 3 - 3 ) 

1=1 

et en faisant la somme sur i des equations ( 13.11 ) , on retrouve la condition 
generale d'equilibre d'un systeme autogravitant isotherme (11.21) 

— > 

V g -P = p m {q) g(q) . 

Nous allons deduire une relation simple entre les densites des differents gaz 
de particules. En utilisant les equations (JH) et ( 13.11 ). on obtient 

— V^ln^) = -^g(q) . (3.4) 



On a done 



1 -* 1 -* 

— Vg-(lnp TO J = — Vg-(lnp m J 



rrii m 3 
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Soit p 0i la densite des particules de masse rrii en un point quelconque, par 
exemple au point q = 0. On a alors la relation suivante entre les densites 
partielles 

fe)* - fe)* 

Nous allons determiner l'equation verifiee par les densites partielles. En 
utilisant l'equation de Poisson du champ gravitationnel 

Vq.g = -47r G p m (q) 

et les equations ( 13.31) et ( 13.41) . on obtient l'equation suivante pour chaque 
densite partielle 

— V| (lnp m J = PrnM ■ (3.6) 

Introduisons le rayon vecteur sans dimension A defini par 



q = a\ , a=\j~r — f t~ (3.7) 

47T G m p 

ou m et p sont des constantes arbitraires ayant respectivement la dimen- 
sion d'une masse et d'une densite de masse. On pose pour la densite des 
particules de masse rrii 

p m = p 0l e Xi(X) , (3.8) 
ce qui impose la condition aux limites pour la fonction Xi 

Xi (S) = . (3.9) 

En utilisant les equations ( 13.61) . ( 13.71) et ( 13.81) et en introduisant les parametres 
sans dimension 

m i Poi /o 1 n\ 

Pi = — , y% = — , (3.10) 

m p 
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on obtient les equations suivantes 

i 

Ces equations determinent les configurations d'equilibre d'un melange de 
gaz autogravitants isothermes. En utilisant les relations ( 13.51) . ( 13.91) et 
3.101 ). on trouve que pour tout i, j 

— > —> 



Hi fJLj 



Ainsi les quantites sont independantes de i. On peut done poser que 



toutes les quantites (1 < i < n) sont egales a une fonction x 



= xW ■ (3-12) 



Les n equations ( 13.111 ) se reduisent alors a une seule equation 



V\ X + e MjX(X) = . (3.13) 

i=i 

Dans le cas de la symetrie spherique, les equations ( 13.111) deviennent 



1 jL a 

/ii A 2 dA 



et 1' equation ( 13.131) devient 



iH x2 S) + ^ e '" xi ^°- (3 - l5) 

Les equations ( 13.141) constituent les equations de la sphere isotherme 
avec plusieurs sortes de particules. La premiere condition initiale est 



X*(A = 0) = 
54 



(3.16) 



Pour que les equations soient regulieres en 0, on impose la deuxieme con- 
dition initiale 



^i(A = 0) = . (3.17) 
dA 

Les equations ( 13.111 ) sont covariantes par la transformation d'echelle 
suivante, a savoir que si les fonctions Xj(A) sont solutions de ces equations, 
alors, etant donne une constante C, les fonctions 

X*$)=Xi(C\) + 21nC (3.18) 

sont aussi solutions de ces equations. Cette symetrie a pour origine l'in- 
variance d'echelle du potentiel de Newton, elle est done valable quelles que 
soient les masses des particules. Cependant la transformation ( 13.181) ne re- 
specte pas la condition aux limites ( 13.91) si bien que l'equation ( 13.131) n'est 
pas invariante par cette transformation. 

Connaissant la solution des equations ( 13.141) . on en deduit les den- 
sites et toutes les grandeurs physiques de la sphere isotherme, comme la 
temperature et les pressions partielles. 

3.1.2 Grandeurs physiques 

Nous allons determiner quelques grandeurs physiques de la sphere isotherme 
en fonction des solutions des equations ( 13.141 ). 

Les parametres rji 

On suppose que le systeme est compose de Ni particules de masse rrii 
et qu'il est contenu dans le volume V. Introduisons les parametres sans 
dimension 

V 3 1 

Pour que Ton puisse considerer le systeme comme autogravitant, il faut que 
l'autogravite et l'agitation thermique jouent toutes les deux un role impor- 
tant dans la physique du systeme ; pour cela, il faut que les parametres 
r]i soient de l'ordre de 1. La limite thermodynamique pertinente est done 
Ni —> oo, V — > oc avec —t fini [T^j. Pour rji — > 0, l'agitation thermique 
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l'emporte largement sur l'autogravite et le systeme se comporte comme un 
melange de gaz parfait. Pour r)i — >• oo, l'autogravite l'emporte largement 
sur l'agitation thermique et le systeme s'effondre sur lui-meme. 

Calculons les parametres rji en fonction des solutions des equations 
( 13.141 ). La conservation de la masse de chaque sorte de particules impose 
la relation suivante 



Ni = [ d 3 q p m {q) . 
Jv 



Plagons nous en symetrie spherique, le volume du systeme etant une sphere 
de rayon Q. En utilisant les equations ( 13.71) . ( 13.81 ). ( 13.101 ) et ( 13.191 ). les 

parametres r\f = rji {^f) 3 definis dans l'introduction par la relation (JHJ) 
verifient la relation suivante 

^ = = M f dA'AVM . (3.20) 

Q T A Jo 

D'apres les equations ( 13.101 ). ( 13.141) et ( 13.201 ). on obtient 

A r)f GmM 1 , , ,^ , o . 

g^ = -g7- = " (3 ' 21) 

On a introduit la masse totale du systeme M = YH=i m i^i- 



Les fonctions fa 

Introduisons les fonctions sans dimension fa = definis dans l'intro- 
duction par la relation @ ou Pi est la pression partielle des particules de 
masse m 8 sur la paroi de la sphere. D'apres les equations (JiJ), ( 13.71 ). ( 13.81 ). 
(ETTH et (E2DD, on a 

D'apres l'equation ( 13.21) . la pression totale du gaz sur la paroi vaut 

p = E T /s • (3 - 23) 

i=l 
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Montrons que les fonctions fa, comme fonctions de Ni, N2,..., A^_i, r\\ 



R 



1 ' 



A^ + i,..., N n , sont solutions d'une equation aux derivees partielles du pre- 
mier ordre. La derivee partielle de fa par rapport a rjf 1 est telle que 

f dfi 

( dfi\ I d\ 



D'apres les equations (EH), (E2DD, JS23J et (E23), il en resulte que 



3/< + Mi J] 



et 




L (3/, - I! 



i]f \d\ ) A 
On en deduit que fa obeit a l'equation suivante 



'if 



V? ( dfi \ 3 ^ + E"=i 77; 



(3.24) 



/< w y 3/, - 1 

Pour j]f = 0, le systeme se comporte comme un melange ideal de gaz 
parfaits. On a done f { (N h N 2 , iV^i, r/f = 0, 7V m , N n ) = 1. 

Plagons nous maintenant a un rayon g inferieur au rayon Q de la paroi. 
D'apres l'equation ( 13.71) . a ce rayon q correspond un rayon reduit X q = X. 
On en deduit d'apres les equations (j3j) et (13.81) que la densite de masse des 
particules de masse rrii et la pression partielle au rayon q sont telles que 



nn Pi(q) _ p mi {q) 




(3.25) 



T P 0l p Ql 

Introduisons le contraste partiel du gaz des particules de masse m 8 qui est 
le quotient de la pression partielle au centre de la sphere et de la pression 
partielle sur la paroi de la sphere 

a = SS = e ^ A) - (3 - 26) 
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Le contraste total qui est le quotient de la pression totale au centre de 
la sphere et de la pression totale sur la paroi de la sphere est, d'apres les 
equations (021) , (E3D|) et (CT5Sj) 



p(Q) ~ Er =1 s e»<») ■ ( ^ 7) 



L'energie 

Pour calculer l'energie, nous allons appliquer le theoreme du viriel ( 11.221) 

£ c + £ = 3 PV 

qui lie l'energie cinetique E c , l'energie totale E et la pression totale de paroi 
P ( 13.231 ). En supposant que toutes les particules sont monoatomiques, on a 
E c =\ Y^i=i N% T. D'apres l'equation ( 13.221) . on obtient le resultat suivant 
pour l'energie 

I = 3 E N > (f* - 1) ■ ( 3 - 28 > 

Nous allons maintenant analyser les instabilites du melange de gaz au- 
togravitants isothermes. 



3.1.3 Stabilite 

Les solutions des equations ( 13.111) sont les configurations d'equilibre 
hydrostatique du melange de gaz autogravitants a l'equilibre thermody- 
namique. lis sont deduits de l'approximation de champ moyen et decrivent 
la phase gazeuse . Comme pour le gaz autogravitant compose par une seule 
sorte de particules, l'etude de la stabilite des configurations d'equilibre du 
melange de gaz autogravitants donne des resultats differents suivant que 
l'on se place dans l'ensemble microcanonique ou dans l'ensemble canon- 
ique car les conditions sur les chaleurs specifiques et les compressibilites 
sont differentes dans ces deux ensembles. 

Nous allons calculer la chaleur specifique a volume constant c v = (ff^)y 
et la chaleur specifique a pression constante c p = (ff')p de la sphere 
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isotherme. Effectuons d'abord le calcul de la chaleur specifique a volume 
constant. En utilisant les equations ( 13.201 ) et ( 13.281 ). on trouve 



n 



U - vi 



R 



i=l 



'dfi 



R 



1 

2 



(3.29) 



En utilisant l'equation ( 13.241 ). on obtient 



i=l 



1 3=1 Mi 



3/i-l 



(3.30) 



Effectuons maintenant le calcul de la chaleur specifique a pression cons- 
tante. En utilisant la formule |27j 



N 

\dVJT 

et les equations ( l3~27ID . (E22D et on trouve 



5> 



/« - ^ 



7? / 5 /i 



i? 



/=1 



U 3 Uf, 



En utilisant l'equation ( 13.241 ). on a 



Cr. 



= -2 ^> + ^JVyj^l 



3/i-l 



i = l 



i=l 



3/i-l 



(3.31) 



(3.32) 



Nous allons calculer maintenant la compressibilite isotherme 
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Kt = —y (fp) T - En utilisant les equations ( 13.201 ) et ( 13.221 ). on trouve 



T 

v KT 



1 



i=l 



f + V 



(3.33) 



En utilisant l'equation ( 13.241 ). on a 



T 



i=i 



3/i-l 



(3.34) 



On obtient la compressibilite adiabatique «5 = — y {jp)$ P ar ia relation 

m 



K S = — KT 

Cy 



On trouve, en utilisant les equations ( 13.291) . ( 13.311 ) et ( 13.331 ). 



T 
V 



«S = 




■i=\ \Ji Hi \grjR 

En utilisant l'equation ( 13.241 ). on obtient 



V 



i=l v 



(3.35) 



+ 



n R Mi E" i —-2 



3/f-l 



. (3.36) 



Dans 1' ensemble canonique, les configurations d'equilibre stable doivent 
avoir une compressibilite isotherme positive et des chaleurs specifiques pos- 
itives. 
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Dans l'ensemble microcanonique, les configurations d'equilibre stable doi- 
vent avoir une compressibilite adiabatique positive. 

Nous avons etudie les configurations d'equilibre de la sphere isotherme 
du melange d' hydrogene et d'helium [15j. II serait interessant de calculer 
les chaleurs specifiques et les compressibilites pour etudier la stabilite de 
ce systeme. 



3.1.4 Lois d'echelle 

Nous avons montre que les melanges de gaz autogravitants isothermes 
constitues de deux sortes de particules obeissent a des lois d'echelle sur la 
masse d'une sphere de rayon q [15j 

M(q) ~ q dH . 

Pour rji = r]2 = (hautes temperatures), on retrouve le gaz parfait (dn = 
3). Pour chaque valeur de les fonctions fi et fi atteignent leurs valeurs 
maximales pour un couple (771,772). Ces couples (r/f, 77^) definissent une 
ligne de point critique ou la chaleur specifique a volume constant diverge. 
On trouve que cLh ~ 1-6 pour tous les points de cette ligne de point critique ; 
cette valeur est done independante de e'est a dire de la composition 
du melange. Ceci manifeste 1' "universalite" des proprietes du systeme a 
l'approche du regime critique. 

Nous avons presente dans cette section l'hydrostatique des melanges 
des gaz autogravitants isothermes. Nous allons etudier dans la prochaine 
section la mecanique statistique des gaz autogravitants composes de deux 
sortes de particules. 



3.2 Mecanique statistique 

Presentons la mecanique statistique des systemes autogravitants com- 
poses de deux sortes de particules dans l'ensemble canonique |T^]. Lorsque 
les nombres Ni et A^2 des deux sortes de particules tendent vers l'infini, 
la fonction de partition est approchee par une integrale fonctionnelle sur 
les deux densites de particules. Le poids statistique des densites est l'ex- 
ponentielle d'une "action effective". On applique l'approximation de point 
col qui est l'approche du champ moyen. Les points col correspondent aux 
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configurations d'equilibre hydrostatique dont les densites de masse par- 
tielles obeissent aux equations ( 13.111) . L'approche du champ moyen decrit 
la phase gazeuse dans la limite thermodynamique autogravitante ou les 
nombres N\ et Ni des deux sortes de particules et le volume V tendent 
vers 1'infini et ou et — f sont finis. Elle montre que, dans cette limite 
thermodynamique diluee, le systeme se comporte comme un melange des 
gaz autogravitants isothermes en equilibre hydrostatique. Elle montre aussi 
que le systeme obeit localement aux equations d'etat des gaz parfaits (Hj). 
Dans l'approche hydrostatique, les equations d'etat ne sont pas determines 
et elles doivent etre supposees. 



3.2.1 Fonction de partition 

On considere un gaz autogravitant dans un volume V et place dans un 
bain thermique a la temperature T constitue de Ni particules de masse mi 
et de Ni particules de masse mi. Les particules de masse m\ et les particules 
de masse mi exercent respectivement sur la paroi une pression partielle P\ 
et une pression partielle Pi. Soit q^\ et pi^ (1 < i < Aq) les positions et 
les impulsions des particules de masse mi et q\^ et p^i (1 < i < Ni) les 
positions et les impulsions des particules de masse mi. Le hamiltonien du 
systeme est 



H = E c + Ep 

Ni -2 N 2 -+i 

^ 2 mi ^2m 2 



G m x 2 G 171 



Ep = - v - v 

_ ^ Gmimi (3 37) 



l<i<N 1 ,l<j<N 2 



W,l ~ q 3 ,2\A 



ou E c est l'energie cinetique et Ep est l'energie potentielle. Notons que 
cut-off defini dans le deuxieme chapitre ( 12.21) et introduit dans Ep permet 
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eviter la divergence de la fonction de partition 



z - wm J l\ w J l\ W e ■ (3 - 38) 

En calculant les integrates gaussiennes sur les impulsions et en introduisant 
les positions sans dimension des particules f^i = et = ^r, on 
trouve que la fonction de partition Z est egale au produit de la fonction 
de partition Zqp du melange ideal de gaz parfaits (GP) de temperature 
T et de volume V contenant Aq particules de masse mi et N2 particules 
de masse mi, par une integrale sur les positions des particules Z int qui 
contient l'information sur l'interaction gravitationnelle : 



Z — ZqP Z int , 

1 / m\ 
Aq~! V~2^r 




/ II d \i d3 r K2 e m Ull+112 ^+V» 12 _ (3.39) 

J Kl<N 1 ,l<k<N 2 

Les parametres rji et 772 sont dermis par les equations 

G m\ Nj G m\ N 2 

m = ^vr ' m = ^w~ (3 - 40) 

et un, U22 et uu sont respectivement l'energie potentielle d'interaction des 
particules de masse m\ entre elles, l'energie potentielle d'interaction des 
particules de masse rri2 entre elles et l'energie potentielle d'interaction des 
particules de masse m\ avec les particules de masse 1112 



u n 
U22 

et 



= - E f , 

_/y n 1 — r,'iL 



= - E 
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Connaissant la fonction de partition, on peut en deduire l'energie libre 
F = — TlnZ et toutes les grandeurs physiques. 



3.2.2 Grandeurs physiques 

L'energie libre F = —TlnZ est, d'apres 1' equation ( 13.391 ) 

F = Fqp — T^ NljN2 ( m , m ) (3.41) 

ou Fqp = — T In Zqp est l'energie libre du melange ideal de gaz parfaits 
de temperature T et de volume V contenant Ni particules de masse m\ et 
A^2 particules de masse JS] • Dans la limite N\ — > oo et N<i — > oo , il est 
bien connu que JHI] 



Fqp = -JViTln 



eV ( m{T\ 



N 2 T In 



eV ( m<{F 



N, 



2 7T 



D'apres les equations ( 13.391) . ( 13.401) et ( 13.411) . la pression 



P 



(8F\ 

\dV )T,N U N 2 



qui s'exerce sur le systeme est 



/ d$ NljN2 



(3.42) 



Introduisons les grandeurs /i et /2 



fi{rji,N 2 ) = 1 - — 



3 \ 9?7i / 

AWlW) - (3.43) 

Dans la limite r}\ — > et 772 — ► (gaz parfait), on a /i(r/i = 0, N2) = 1 et 
f2(Ni, 7)2 = 0) = 1. En integrant la relation ( 13.431) . on obtient 
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= 3 JVi /"" dx 1 ^ ^ 

+ 3 iV 2 " f dz izWl . (3.44) 
Jo x 

On peut exprimer toutes les grandeurs thermodynamiques en fonction 
des grandeurs f\ et / 2 grace a cette equation. D'apres l'equation ( 13.411 ). 
l'energie libre est 



F 



F GP - 3 N lT £ dx l^lM 
Jo 



(3.45) 



L'energie moyenne est en utilisant la relation < £" >= F — T (§f?)y 



< E > = 3 N\T 



h 



i 



+ 3iV 2 T 



/2 



1 



On en deduit que l'entropie S = s'exprime suivant 



(3.46) 



S = S GP + 3 N x 
+ 3iV 2 



Jo 



/2 _ 1+ r ^-f^ 

Jo 



X 



OU 



(3.47) 



\ v 3 i 
In 1 — m 

M 2 



/ miT\ 
V 2tt J 



5' 

+ 2 



+ iV 2 



V 3, /m 2 T\ 5 
n — + - In + - 

N 2 2 V 2tt y 2 



est l'entropie du melange ideal de gaz parfaits de temperature T et de 
volume V contenant N\ particules de masse mi et iV 2 particules de masse 
m 2 [31]. 

Nous allons maintenant exposer l'approximation de champ moyen qui, 
nous allons le voir, conduit a l'hydrostatique que nous avons presentee dans 
la premiere partie de ce chapitre et qui est exacte dans la limite N± — ► oc 
et N 2 — ► oc. 
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3.2.3 Champ moyen 



Dans la limite thermodynamique autogravitante (N\ oo, N 2 ^ oc 
et V — > 00 avec — r et — | finis), on applique l'approximation de champ 
moyen. L'integrale Zi nt devient une integrale fonctionnelle sur la densite 
Pi{r) des particules de masse mi et la densite P 2 {t) des particules de masse 
m 2 



Zint = f D Pl (.) D P2 {.) ^ ^ exp[-s c ( Pl (.),p 2 (.),b h b 2 )} (3.48) 



avec 1' "action effective" 



771 /" d 3 r d 3 f / , 
- iVi — / -p — ^- pi(r) pi(r ) 
z J |r — r I 

7/ 2 /" d 3 f d 3 f' , 
N 2 — / — — — p 2 (r) p 2 (r ) 
z J \r — r I 

T=J — =rrr £>i(r) /o 2 (r 
|r — r I 

iVi / d 3 r pi(r) ln/?i(r) 
7V 2 y d 3 f /9 2 (r) ln/? 2 (r) 



y d 3 r p\{f] 



ib 2 N 2 



d 3 r p 2 (r) 



(3.49) 



Les integrations sur la position r se font sur le volume unite. L'integrale 
fonctionnelle ( 13.481) est dominee pour Aq — > 00 et N 2 — > 00 par le point col 
de 1' " action effective" ( 13.491) qui verifie les relations suivantes 



(Pcol,h Pool, 2, bcol,l, bcol,2) = , Trr(Pcol,l, Pcol,2, b co l,l, b C ol,2) = , 



dbi 



db, 



Ss c &s c 

(Pcol,h Pcol,2, b co i 7 i, b co i j2 ) = , - - - ( ^ {Pcol,l, Pcol,2, bco^i, b co i^ 2 ) = 



Spl(-) 



SP2(.) 
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Les deux premieres relations imposent la normalisation des deux densites 

J dV p coll {r) = 1 , J d¥ Pcol , 2 {r) = 1 . (3.50) 

Les deux dernieres relations imposent que les deux densites soient solutions 
des equations de point col suivantes 



-rz — —\ Pcoi,i(r) 
\r — r | 

-rz — —\ Pcoi,i(r) + a coZj i 
\r — r | 

inPcoi,2{r) = - m — —\ Pcoi,i{r ) 
fi J \r — r | 

— Pooler ) + a co i }2 ■ (3.51) 
\r — r | 

a C oi,i = i b co ^i — 1 et a co i^ = ib co i,2 — 1 sont les multiplicateurs de Lagrange 
associes a la condition de normalisation des deux densites ( 13.501) . On a 
introduit le rapport fi = En appliquant le Laplacien aux equations de 
point col et en introduisant les fonctions $i(r) = ln/? co ^i(f) et <3>2(r) = 
In Pooler), on trouve 



V£ $i(f) + 4 7T 7/i e* 1 ^ + 4 7T /i r/2 e*™ = 

V£ $ 2 (f) + 4 7T i 7/! e* 1 ^ + 4 7T 7/2 e* 2 ^ = . (3.52) 
fi 

Les densites sans dimension p co i,i{r) = e* 1 ^ et p co i,2{r) = e* 2 ^ sont liees 
aux densites de masse p mi {q) et p m2 (q) introduites dans la premiere partie 
de ce chapitre par les relations 



miNi mi Ni $ m 

Pm 1 {q) = PcolAn = —y- e > 

77^2 7772JV2 ^ 

f J m 2 \Hj y Pcol,2\i ) y e J 

(f = Wf. (3.53) 
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En utilisant les equations ( 13.401) et ( 13.531 ). on trouve que les equations de 
point col sont identiques aux equations d'equilibre hydrostatique ( 13.111) 

^-V|lnp mi + \pmM + PmM\ = 

et 

— V|lnp m2 + [p mi (q) + PmM] = 0. 

777-2 i 

Les solutions de point col sont done identiques aux configurations d'equi- 
libre hydrostatique des melanges des gaz autogravitants isothermes. Dans 
la limite thermodynamique autogravitante (Aq — ► oo, A^ — > oo et V — > 

oo avec — r et — f finis), la mecanique statistique montre que le systeme 

\/3 y 3 



obeit aux equations d'etat des gaz parfaits (J2J) et aux equations d'equilibre 
hydrostatique ( 13.111) . 

D'apres l'equation ( 13.531 ). la pression partielle des particules de masse 
mi et la pression partielle des particules de masse au point r sont 

PM = Pcoi,i(rl , P 2 (q) = Pco^r) ■ (3-54) 
Dans le cas de la symetrie spherique, les fonctions fi = et f 2 = j^f 



(eqs. ( 13.431 )) coincident avec les fonctions f\ et / 2 , introduites dans les 
equations ( 13.221) dans le cadre de l'hydrostatique. Rappelons que ces fonc- 
tions fi et f 2 obeissent aux equations aux derivees partielles ( 13.241) . On 

a 

PiV = f P2V = , 
NiT h ' N 2 T h ' 

Les equations ( 13.241) s'integrent de cette maniere 

3 / — [l-/i(^^2)] = 3[/i(^,^)-l]+(^ + ^)-ln/ 1 ( ?7 f,^ 
J x 

et 



3 f 2 — [1 - f 2 (N 1 , x)} = 3 [f 2 (N h x) - 1] + ( ^ + 7/) - ln/ 2 (Aq, r?f 
Jo x \ A* / 



(38 



En utilisant les equations ( 13.411) et ( 13.471) . on en deduit que l'energie libre 
et l'entropie verifient les relations suivantes 



et 



= ^i[3(l-/i)-(r/f + M^) + ln/i] 



+N 2 



3(1-/2)- ^ + ^)+ln/2 



(3.55) 



S-S GP = Aq[6 (/i-l) + (r/f + rf)-Wi] 



+N 2 

En generalisant les relations 
on a 



R 



6(/ 2 -l)+f^ + ^ 



ln/2 



(3.56) 



et ( 13.561) pour n sortes de particules, 



F — Fqp 
T 



i=l 



et 



R 



3 (1 - fi 



(3.57) 



6(/i-l)+A*<X) — 
i=i ^' 

II serait interessant d'etudier la validite du champ moyen, en calculant 
les petites fluctuations autour du point col comme ceci a ete fait dans le cas 
avec des particules identiques [7j. Des calculs Monte Carlo nous permet- 
traient de verifier que le champ moyen decrit bien la phase gazeuse ; ils nous 
permettraient aussi de determiner la transition de la phase gazeuse vers 
la phase collapsee. On pourrait comparer les resultats des calculs Monte 
Carlo avec les resultats sur la stabilite des configurations d'equilibre hy- 
drostatique deduits dans la premiere partie de ce chapitre a partir du com- 
portement des chaleurs specifiques et des compressibilites, pour voir s'ils 
coincident. 

II serait aussi interessant d'etudier la mecanique statistique des systemes 
autogravitants composes de plusieurs sortes de particules dans l'ensemble 
microcanonique et dans l'ensemble grand-canonique. 



R 



(3.58) 
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Chapitre 4 



Systemes autogravitants en presence 
de la constante cosmologique 



De recentes observations astrophysiques ont mis en evidence que l'u- 
nivers est rempli d'une energie noire que Ton modelise par la constante 
cosmologique A des equations d'Einstein de la relativite generale. Nous 
presentons la limite non relativiste des equations d'Einstein avec la cons- 
tante cosmologique et nous etudions la mecanique statistique des systemes 
autogravitants en presence de la constante cosmologique. Pour un systeme 
autogravitant en presence de la constante cosmologique A, comportant 
A" particules dans un volume V, l'energie thermique est de l'ordre de N, 
1' energie autogravitante (due aux seules interactions gravitationnelles entre 
particules) est de l'ordre de et l'energie de la constante cosmologique est 
de l'ordre de AV. Pour que ces trois energies soient de meme importance, 

1 2 

il faut que A" ~ V s et A ~ V~z. La limite thermodynamique pertinente 
est done Af^oo,y^ocetA^0 avec -^r et AVz finis. Dans cette lim- 
ite thermodynamique, nous avons montre que l'approche du champ moyen 
decrit exactement la phase gazeuse [19j. Nous avons faits des calculs Monte 
Carlo dans l'ensemble canonique [20] ; ils confirment que le champ moyen 
decrit tres bien la phase gazeuse et qu'il cesse d'etre valable lorsque la com- 
pressibilite isotherme diverge en devenant negative ; une transition s'opere 
alors de la phase gazeuse vers la phase collapsee. Les calculs Monte Carlo 
permettent d'etudier la phase collapsee. 



71 



4.1 La constante cosmologique 



Nous presentons dans cette section la constante cosmologique et la limite 
non relativiste des equations d'Einstein avec la constante cosmologique. 

4.1.1 La constante cosmologique et l'energie noire 

La constante cosmologique a ete introduite par Einstein dans les equa- 
tions de la relativite generale en 1917 [27]. Elle permet d'obtenir un univers 
homogene statique comme solution de ces equations. Les observations du 
redshift de la lumiere emise par les galaxies ont mis en evidence que les 
galaxies s'eloignent les unes des autres et que l'univers est actuellement en 
expansion [38]. Des solutions des equations de la relativite generale sans 
la constante cosmologique decrivent un univers homogene dynamique en 
conformite avec ces observations jSSJ- Cependant, de recentes observations 
basees sur la mesure de la lumiere d'un type de supernovae de luminosite 
intrinseque pratiquement uniforme ont mis en evidence qu'aujourd'hui les 
trois-quart de l'energie de l'univers ne sont pas constitues de matiere et 
de rayonnement mais d'une energie appelee energie noire (TBI El HBJ- 
Ses proprietes observees actuellement sont proches de celles modelisees 
par la constante cosmologique. La constante cosmologique A agit comme 
une densite d'energie uniforme extremement faible, ayant un effet repulsif 
sur la matiere ; elle a pour effet d'accelerer l'expansion de l'univers. Son 
importance vient du fait qu'elle remplit tout l'univers. 

4.1.2 Gravitation non relativiste 

Presentons les modifications qu'apporte la constante cosmologique a la 
gravitation non relativiste. Les equations de la relativite generale sont jlDJ 



ou est le tenseur de Ricci et T@ est le tenseur d'energie-impulsion. En 
considerant que la matiere est non relativiste, c'est a dire que sa pression est 
negligeable devant sa densite de masse p m , le tenseur d'energie-impulsion 
en presence de la constante cosmologique A a pour expression 




(4.1) 



72 



T* = Pm S a0 5 m + A € ■ 

Dans la limite non relativiste, la composante 00 des equations de la rela- 
tivite generale ( 14.11 ) devient [19j 

Vq.g + 4ttG (p m (q) - 2A) = . (4.2) 

II s'agit de l'equation reliant au point de position (fie champ gravitationnel 
g a la densite de masse p m en presence de la constante cosmologique A. Pour 
A = 0, on retrouve l'equation de Poisson de la gravitation newtonienne 
( 11.41 ). Le champ gravitationnel s'integre de la maniere suivante 

g {ql = -G I d 3 q'p m (t) ^4 + ^ 9 , (4.3) 

J \q-q r 3 

le premier terme attractif etant la contribution de la matiere et le second 
terme repulsif etant la contribution de la constante cosmologique. La cons- 
tante d'integration est choisie nulle dans l'equation ( 14.31) . ce qui revient 
a eliminer les champs gravitationnels exterieurs en placant le centre de 
masse au point q = 0. Le potentiel gravitationnel au point q defini par 
g(q) = — VqV vaut 

d 3 q'p m (q) _ 4ttGA_ 2 

\q-q\ 3 

Le hamiltonien d'un systeme autogravitant en presence de la constante 
cosmologique dont les particules ont comme masses mi, m^v, comme 
positions q\ : ... : qN et comme impulsions p\, est 




H = E c -\- Ep 

Ep - -G ^ „ry - — ^ m q, , (4.4) 

l<i<j<N ^ ^ n i=l 

E c etant l'energie cinetique et Ep l'energie potentielle. 
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Presentons maintenant l'hydrostatique des gaz autogravitants en pre- 
sence de la constante cosmologique qui est deduite de 1' approximation de 
champ moyen de la mecanique statistique. 

4.2 Hydrostatique 

Le systeme autogravitant en presence de la constante cosmologique ex- 
iste sous deux phases qui ne peuvent pas coexister ensemble, une phase 
gazeuse et une phase collapsee. Nous allons presenter l'hydrostatique d'un 
gaz autogravitant isotherme en presence de la constante cosmologique qui 
est deduite de l'approche du champ moyen de la mecanique statistique. 
Dans la limite thermodynamique ou le nombre de particules N, le volume 
V et la constante cosmologique A verifient iV^oo,"l/^ocetA^0 
avec pr fini et AV^ fini, l'approche du champ moyen decrit exactement la 
phase gazeuse du systeme. 

4.2.1 Equilibre hydrostatique 

Rappelons qu'un gaz autogravitant est en equilibre hydrostatique si les 
forces de pression et les forces gravitationnelles se compensent, ce qui se 
traduit par la relation ( 11.21 ) 

-V g -P + p m (q) g(q) = 

entre la pression P du gaz au point q, la densite de masse p m au point q 
et le champ gravitationnel g engendre par le gaz autogravitant au point q. 
Son expression est donnee par l'equation ( 14.31) lorsqu'il est en presence de 
la constante cosmologique. En utilisant l'equation du champ gravitationnel 
( 14.21 ). on obtient la relation suivante qui est la condition d'equilibre hydro- 
statique d'un gaz autogravitant en presence de la constante cosmologique 

V«f (— VfP) = -4ttG (p m (q) - 2A) . (4.5) 

4.2.2 Equilibre thermodynamique 

Dans la limite thermodynamique (N ^oo,y^ooetlaA^0 avec 
et AV 3 fini), l'approche du champ moyen de la mecanique statistique 



74 




0.2 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 ^ 3 

Fig. 4.1 - La courbe f(rj,R\) pour R\ = 0,0.3,1,1.5 par l'approche du champ moyen 

ou / = (voir equation lJOT|l ), rj = G jf* = rj R (J^) 5 (voir equation J32D ) 

et i?A = (voir equation 0). Les configurations d'equilibre stable dans l'ensem- 
ble canonique sont representees par les points compris entre le point (j] — 0, / = 1) 
et le point (rj = rj (R A ) , f (R A ) = f \r] (R A ) , Ra)) represents par un +. Pour ces con- 
figurations, les chaleurs specifiques et la compressibility isotherme sont positives. Au 
point (i] (R A ), f (R\)), la chaleur specifique a pression constante et la compressibilite 
isotherme divergent en devenant negatives ; il s'agit du point d'instabilite de la sphere 
isotherme dans l'ensemble canonique. On a rj (R A = 0) = 1.510..., t] (Ra = 0.3) = 1.63..., 
r) (R A = 1) = 2.04... et r] (R A = 1.5) = 2.55... . 
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decrit exactement la phase gazeuse. Elle montre que localement la pression 
et la densite de masse obeissent a l'equation d'etat des gaz parfaits (JH) 

T 

p(q) = —pm{q) 

m 

ou m est la masse de chacune des particules du gaz qui sont supposees ici 
identiques et ou T est la temperature constante. A partir de cette equation 
et de l'equation ( 14.51) . on obtient l'equation de la densite du gaz autograv- 
itant a l'equilibre thermodynamique 

V|(lnp TO ) + (p m (q) - 2A) = . (4.6) 

En symetrie spherique, les grandeurs physiques ne dependent que de la 
distance q par rapport au centre de la sphere. L'equation ( 14.61 ) devient 

~ 2 ^ {q 2 lnp m ) + (p m (q) - 2A) = . (4.7) 

II s'agit de l'equation de la sphere isotherme en presence de la constante 
cosmologique. La constante cosmologique a brise la covariance de l'equation 
de la densite par la transformation d'echelle ( 11.151) . On ne peut plus reduire 
cette equation du second ordre a une equation du premier ordre comme 
c'est le cas lorsque A = 0. 



4.2.3 Les parametres 77 et £ 

Considerons un systeme autogravitant isotherme de temperature T en 
presence de la constante cosmologique A. II est compose de particules 
de masse m dont la masse totale M = mN est contenue dans un volume 
V. Ce systeme a deux parametres significatifs sans dimension. Le premier 
parametre est 

GmM Gm 2 N 

rj = 1 = 1 4.8 

yi T T 

defini dans le chapitre 1. On rappelle que le parametre ij est le quo- 
tient de deux energies caracteristiques d'une particule en interaction avec 
l'ensemble du systeme autogravitant. Ces deux energies caracteristiques 
sont i qui est de l'ordre de son energie gravitationnelle d'interaction 

V3 00 
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avec 1' ensemble du systeme autogravitant et T qui est de l'ordre de son 
energie cinetique. Le deuxieme parametre est 



Z = f • (4-9) 

Ce parametre est le quotient entre deux energies caracteristiques d'une 
particule en interaction avec la constante cosmologique. Ces deux energies 

2 

caracteristiques sont 2G m AVs qui est de l'ordre de son energie grav- 
itationnelle d'interaction avec A et T qui est de l'ordre de son energie 
cinetique. 

Le rapport entre ces deux parametres 

est le rapport deja defini dans l'introduction (eq.(J7j)) entre la quantite 
d'energie de la constante cosmologique et la masse de matiere dans le 
systeme ; il exprime l'importance relative de la constante cosmologique et 
de la matiere dans le systeme. 

La limite thermodynamique pertinente des systemes autogravitants en 
presence de la constante cosmologique AestAf^oc,V^oc,A^O 
avec -jjr et AVs finis. Pour cette limite, 77 et £ sont d'ordre unite (Ra ~ 1) 
et l'autogravite des particules (les interactions gravitationnelles mutuelles 
des particules ) et la constante cosmologique jouent toutes les deux un role 
important. Lorsque 77 et £ tendent vers 0, l'energie cinetique du systeme 
l'emporte largement sur son energie gravitationnelle et le systeme se com- 
porte comme un gaz parfait. Lorsque £ — > et 77 est fini (it! a ~ 0), les effets 
de la constante cosmologique sont negligeables devant ceux de l'autogravite 
des particules et on retrouve les systemes autogravitants etudies dans les 
premiers chapitres. Le cas ou les effets de l'autogravite des particules sont 
negligeables devant ceux de la constante cosmologique (77 — > et £ est fini, 
R\ ^> 1) est presente plus loin dans ce chapitre. 

Exprimons la difference 77 — £ a partir des solutions de l'equation ( 14.61) 
dans le cas de la sphere isotherme. A partir de l'expression de la masse 
contenue a l'interieur du volume V 



M = m N = / d 3 q p m (q 
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de 1' equation ( 14.61 ) et du theoreme de Green-Ostrogradski sur le volume V, 
on trouve que 

/-» 4:7tG m 

dSVq lnp m + {mN -2AV) = 0, 

la surface d'integration etant la paroi entourant le systeme. Dans le cas de 
la symetrie spherique ou le volume est une sphere de rayon Q, la difference 
des parametres T] R = 77 (4^) 3 et £ R = £ (^) 3 est 

rf-t* = -Q(±Qn*.))^. (4.10) 

La presence de la constante cosmologique a pour effet de transformer 
l'equation ( 11.161) en l'equation ( 14.101) en substituant le parametre r\ , pro- 
portionnelle a la masse de matiere M, par la difference r] R — £ , propor- 
tionnelle a la difference entre la masse de matiere M et l'energie 2AV de 
la constante cosmologique. 



4.2.4 Densite de la sphere isotherme 

Nous allons considerer un gaz autogravitant a l'equilibre thermody- 
namique en symetrie spherique (sphere isotherme). Le systeme est contenu 
dans une sphere de rayon Q sur la paroi de laquelle est exercee la pression 
P. Nous allons etudier les variations de la densite de massse en fonction 
du rayon q (0 < q < Q) [T$] . En l'absence de la constante cosmologique, 
la densite de matiere de la sphere isotherme est toujours une fonction 
decroissante du rayon q, ceci est du a 1' effet attractif de l'autogravite des 
particules. En presence de la constante cosmologique, il y a competition 
entre 1' effet attractif de l'autogravite des particules et 1 'effet repulsif de 
la constante cosmologique. Pour les configurations d'equilibre stable de la 
sphere isotherme, les variations de la densite ont trois comportements : 

-pour i?A < 1, la densite est une fonction decroissante du rayon q. Les 
effets de l'autogravite l'emportent sur ceux de la constante cosmologique. 

-pour i?A = 1, la densite est une fonction uniforme du rayon q. Les 
effets de l'autogravite et de la constante cosmologique se compensent et 
le systeme se comporte comme un gaz parfait. Le systeme decrit l'univers 
homogene et statique d'Einstein dans une version non relativiste. 
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-pour R A > 1, la densite est une fonction croissante du rayon q. Les 
effets de la constante cosmologique l'emportent sur ceux de l'autogravite. 
Introduisons la grandeur sans dimension / (fig J4.ll) 

MRa) = = ^A»W) (4-11) 

en rappelant que la pression sur la paroi P et la densite sur la paroi p m (Q) 
sont liees localement par l'equation d'etat des gaz parfaits (DP). 



4.2.5 Stabilite de la sphere isotherme 

Nous avons determine la stabilite des configurations d'equilibre de la 
sphere isotherme, solutions de l'equation ( 14. 71) . dans l'ensemble canonique 
[20J. Etant donne une valeur de R A , ces configurations d'equilibre sta- 
bles correspondent aux points du diagramme de phase (fig J4.ll) du point 
(?7 = 0,/ = 1) au point (rj = rj (R A ) , f (R A ) = f{rj {R A ), R A )). Pour 
ces configurations, les chaleurs specifiques et la compressibilite isotherme 
sont positives. Au point (r] (R A ), f (R A )), la chaleur specifique a pression 
constante et la compressibilite isotherme divergent en devenant negatives. 
Cette configuration est le point d'instabilite de la sphere isotherme dans 
l'ensemble canonique. Le parametre t] est une fonction croissante de R A , 
la presence de la constante cosmologique a pour effet d'etendre la sta- 
bilite du gaz autogravitant. Par exemple, on a r] (R A = 0) = 1.510..., 
rj (R A = 0.3) = 1.63..., r) (R A = 1) = 2.04... et rj (R A = 1.5) = 2.55... . 

Presentons maintenant la mecanique statistique des systemes autograv- 
itants en presence de la constante cosmologique. 



4.3 Mecanique statistique 

Nous avons developpe la mecanique statistique des systemes autogravi- 
tants en presence de la constante cosmologique dans l'ensemble canonique. 
L'approche du champ moyen decrit exactement la phase gazeuse dans 
la limite thermodynamique autogravitante (N ^oo,y^ooA^0 
avec — r et AV^ finis). Le champ moyen montre que le systeme obeit a 
l'equation ( 14.61) de l'hydrostatique et permet d'obtenir l'equation d'etat 
du systeme, celle-ci correspondant localement a l'equation d'etat des gaz 
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parfaits ((H) [19|. En revanche, en hydrostatique l'equation d'etat n'est pas 
derivee et doit etre supposee. 



4.3.1 Fonction de partition 

Considerons un gaz de N particules de masse m dans un volume V. II est 
place dans un thermostat a la temperature T et une pression P s'applique 
sur la paroi enfermant le systeme. Les particules interagissent entre elles 
par la gravite et sont en presence de la constante cosmologique. La fonction 
de partition est 



m J l\ (2^)3 



avec 



H = E c + U , 



N ^ 2 



^ 2m 



i=l 



£ - = _ £ ^7 - -3- £ a 2 • ( 4 - 12 > 

l<i<j<N |y yjl i=l 



etant le hamiltonien de l'equation ( 14.41) avec des particules identiques. En 
calculant les integrates gaussiennes sur les impulsions et en introduisant les 
positions sans dimension des particules f\ = on trouve que la fonction 
de partition Z est egale au produit de la fonction de partition Zqp du gaz 
parfait (GP) de temperature T et de volume V contenant N particules de 
masse m par une integrale sur les positions des particules Z int qui contient 
l'information sur l'interaction gravitationnelle et la constante cosmologique 




'volume unite 



Y\ d 3 n e v "p^'-'^'+f £ u N (f 1: ...,r N ) _ (4.13) 
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Les parametres 77 et £ sont dermis par les equations ( 14.81 ) et ( 14.91) . up est 
l'energie potentielle sans dimension due a l'autogravite des particules et un 
est l'energie potentielle sans dimension due a la constante cosmologique 



u P (f h ...,r N ) = i ^ 1 , ujv(fi r .,fjv) = X)r 



TV 

2 



7V ^ 7*' — 7*' 

l<i<j<N 1 1 J|a i=l 



4.3.2 Champ moyen 

Dans 1' approximation de champ moyen, l'integrale ( 14.131 ) devient 
une integrale fonctionnelle sur la densite p(r) 

^<n* = J Dp(.) ^ exp[-iVs c (p(.),&)] (4.14) 
avec 1' "action effective" 



s c(p(-),b) = J d 3 fp(f) lnp(r) 

— <$; / d 3 r p(f) r 2 + 26 



: -4 



z / r — r 



1 



d 3 r p(r) 



.(4.15) 



L'integrale fonctionnnelle ( 14.141) est dominee pour N — ► 00 par le point 
col de 1' "action effective" ( 14.151) qui verifie les relations suivantes 

-^{Pcoh bed) = , JjTjO>ccl, hoi) = . 

La premiere relation impose la normalisation de la densite 



J d 3 f Pcol (r) = 1 



(4.16) 



La seconde relation impose que la densite soit solution de l'equation de 
point col 



In p co i(r) - 77 J 



d 3 r 



PcoiKT 



2^r , - 2 



(4.17) 



&coi — ib co i — 1 etant un multiplicateur de Lagrange associe a la condition de 
normalisation de la densite ( 14.161 ). En appliquant le Laplacien a l'equation 
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de point col ( 14.171 ) et en introduisant la fonction $(r) = In p co i(r), on trouve 
que 

V|$(r) + 4 7r(7/e $ ^ - A =0. (4.18) 

La densite sans dimention p co i(f) = e $ ( r ) est liee a la densite de masse 
p m (q) introduite en hydrostatique par la relation 

TYl N l 

Pm(q) = -y- Pcdif) , ?=V*f. (4.19) 

En utilisant les equations ( 14.81 ). ( 14.91) et ( 14.191 ). on trouve que l'equation de 
point col est identique a l'equation d'equilibre hydrostatique ( 14.61) . Les so- 
lutions de point col du systeme sont done equivalentes aux configurations 
d'equilibre hydrostatique du systeme autogravitant isotherme avec cons- 
tante cosmologique. Dans la limite thermodynamique (N — > oo, V —>■ oo 
et A — > avec et AV 3 finis) , la mecanique statistique montre que le 
systeme obeit a l'equation d'equilibre hydrostatique ( 14.61) et a l'equation 
d'etat locale des gaz parfaits inhomogenes (CD). 

4.3.3 Calculs Monte Carlo 

Nous avons effectue des calculs Monte Carlo (500 < N < 2000) [20J qui 
confirment que la phase gazeuse est bien decrite par le champ moyen et qui 
permettent d'etudier la phase collapsee et la transition de la phase gazeuse 
vers la phase collapsee. Conformement a nos previsions, celle-ci a bien lieu 
lorsque la compressibilite isotherme diverge en devenant negative. Nous 
avons effectue les calculs Monte Carlo dans un cube et dans une sphere, 
ce qui nous a permis d'etudier l'influence de la forme de la paroi sur les 
systemes autogravitants. En presence de la constante cosmologique, les 
resultats sont dependants de la forme de la paroi, alors qu'en l'absence de 
la constante cosmologique, les resultats sont moins dependants de la forme 
de celle-ci. L'action de la constante cosmologique etant plus importante au 
niveau de la paroi qu'au centre du systeme, la forme de celle-ci a d'avantage 
d'influence en presence de la constante cosmologique. 

Presentons maintenant la limite Ra ^> 1 (eq.flZj)) ou la constante cos- 
mologique domine l'autogravite. 
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4.4 Limite R A > 1 



Nous considerons le cas ou la constante cosmologique domine l'auto- 
gravite (Ra ^> 1) Le hamiltonien ( 14.121 ) ou l'energie potentielle auto- 
gravitante a ete negligee devient 

tt Pi 4irGmA ^ ^ 2 
# = 2^2^ • 

i=l i=l 

En symetrie spherique, la fonction de partition s'exprime (sans approxima- 
tion de champ moyen) suivant 

3 N 

Z= V W (^Y e-^» (4,0) 
avec ^ = £ (eq.63) et 

e-< fl ) = 3 j\y y 2 exp (^y 2 ) . (4.21) 
La grandeur / definie par l'equation ( 14.111) vaut 

= e < fl ) e ^ . (4.22) 

D'apres l'equation ( 14.211 ). elle verifie l'equation differentielle du premier 
ordre 

J^-2 (1 - f) + Y (4 ' 23) 

avec la condition limite 



m R = o) = i . 

En utilisant les equations ( 14.201 ) et ( 14.221 ). la fonction de partition s'exprime 
suivant 

M , 3 jV 

N TA 2 t R 



"iVT V~27T 



Z = ^— [ — ] e N ^ f~ N . (4.24) 
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On en deduit l'energie libre F, l'energie E, la pression sur la paroi P, 
l'entropie S, la chaleur specifique c v et la compressibilite isotherme kt en 
fonction de la grandeur / en utilisant les equations ( 14.221) et( !4.24l) 



PV 
NT 



S — Sgp 3 , 
N = 



[l-/(ai-ln/(^) + 



R 



R 



i + 5 /(C 



1? 



et 



-i 



3 



Les grandeurs Fqp et Sgp sont respectivement l'energie libre et l'entropie 
du gaz parfait d'energie E et de volume V compose de particules. Toutes 
les configurations ont une compressibilite isotherme positive, elles sont sta- 
bles dans 1' ensemble canonique. La densite de masse p m s'exprime en fonc- 
tion de la distance q par rapport au centre de la sphere (0 < q < Q) et de 
la constante cosmologique suivant 



Vpm(q) 

mN 



= e 



exp 




A cause de l'effet repulsif de la constante cosmologique, la densite de masse 
est une fonction croissante du rayon q. 

Discutons maintenant de l'importance de la constante cosmologique 
dans les objets astrophysiques autogravitants. 



4.5 Discussions 

L'importance de la constante cosmologique est mesuree par le rapport 
(jZj) entre l'energie de la constante cosmologique et la masse de la matiere. 
Plus un systeme est dense et plus l'importance de la constante cosmologique 
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est faible. A l'echelle de l'univers qui est homogene, ce rapport vaut 4. 
A partir d'une echelle de distance inferieure a lOMpc, l'univers a une 
structure fragmentee en elements de plus en plus denses en matiere. Plus 
la structure etudiee est petite et plus l'importance de la constante cos- 
mologique (qui est repartie de maniere uniforme) est faible relativement a 
la matiere. Elle est negligeable pour le milieu interstellaire mais doit etre 
prise en compte pour les structures de galaxies les plus grandes que sont 
les amas et les superamas. 
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Conclusions et perspectives 



Dans cette these, nous avons etudie la mecanique statistique des sys- 
temes autogravitants comportant plusieurs sortes de particules et la me- 
canique statistique des systemes autogravitants en presence de la cons- 
tante cosmologique. Ces deux contributions sont entierement nouvelles. 
Pour ces deux types de systemes autogravitants, nous avons developpe 
l'approche du champ moyen qui decrit exactement les phases gazeuses de 
ces systemes dans leurs limites thermodynamiques pertinentes respectives. 
Pour les systemes autogravitants comportant plusieurs sortes de particules, 
cette limite thermodynamique est la limite ou les nombres de particules N,, 



i 

Ni 



et le volume V tendent vers l'infini et ou les rapports — f sont finis. Pour 
les systemes autogravitants en presence de la constante cosmologique, cette 
limite thermodynamique est la limite ou le nombre de particules N et le 
volume V tendent vers l'infini ou la constante cosmologique A tend vers 
et ou les rapports et AVs sont finis. Dans leur limite thermodynamique 
respective, l'approche du champ moyen montre que ces systemes autograv- 
itants obeissent aux equations de l'hydrostatique et a une equation d'etat 
qui localement est l'equation d'etat des gaz parfaits. Nous avons calcule 
les grandeurs thermodynamiques de ces systemes. Nous avons analyse leur 
stabilite. Nous avons effectue des calculs Monte Carlo pour le systeme auto- 
gravitant en presence de la constante cosmologique dans l'ensemble canon- 
ique. lis montrent que le champ moyen decrit tres bien la phase gazeuse 
et que la transition de la phase gazeuse vers la phase collapsee a lieu dans 
l'ensemble canonique lorsque la compressibilite isotherme diverge. II serait 
interessant de faire des calculs Monte Carlo dans l'ensemble microcanon- 
ique pour verifier, conformement a nos previsions que la transition de phase 
arrive dans cet ensemble lorsque la compressibilite adiabatique diverge. II 
faudrait egalement faire des calculs Monte Carlo pour les systemes auto- 
gravitants comportant plusieurs sortes de particules. Nous avons montre 
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que les systemes autogravitants comportant plusieurs sortes de particules 
obeissent a des lois d'echelle sur leur masse. Le milieu interstellaire qui 
est compose de plusieurs sortes d'atomes et de molecules a done ses lois 
d'echelle sur sa masse reproduites par les systemes autogravitants compor- 
tant plusieurs sortes de particules. Par son action repulsive, la constante 
cosmologique augmente la stabilite des systemes autogravitants. 

Nous avons etudie dans cette these la mecanique statistique et l'hydro- 
statique des systemes autogravitants. II serait interessant de developper 
l'hydrodynamique des systemes autogravitants qui verifient les equations 
(IB. II ) et ( IB. 21) . Dans l'annexe B, nous exposons la theorie des perturbations 
par rapport a un fluide autogravitant homogene statique. II serait utile 
d'etudier la theorie des perturbations avec comme ordre zero un gaz au- 
togravitant en equilibre hydrostatique. Ce serait un autre moyen d'etudier 
la stabilite des gaz autogravitants en equilibre hydrostatique ; ce serait 
aussi un moyen d'etudier leur collapse. L'hydrodynamique des systemes 
autogravitants est actuellement etudiee en cosmologie pour montrer com- 
ment se forment les structures a partir d'un fond homogene de matiere 
dans l'univers en expansion. La theorie des perturbations par rapport au 
modele d'Einstein-de Sitter J2BJ qui est le modele d'univers plat en expan- 
sion domine par la matiere a ete etudiee jH], S2j • Pour les temps infinis, 
cette theorie des perturbations diverge, il serait interessant d'appliquer le 
groupe de renormalisation dynamique |43j pour explorer le comportement 
de la theorie pour les temps infinis et prevoir ainsi revolution des structures 
dans l'univers. 
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Annexe A 



Gaz autogravitants polytropiques 

Nous allons presenter brievement les gaz autogravitants polytropiques 
qui jouent un role important dans la comprehension de la physique des 
etoiles pEJ Ej ESI EH HHl US]- Nous allons tout d'abord presenter les trans- 
formations polytropiques que subissent ces systemes. 

A.l Les transformations polytropiques 

Une transformation polytropique jlj est un transformation thermody- 
namique ou la variation de chaleur est proportionnelle a la variation de 
temperature. En considerant une transformation polytropique infinitesi- 
male, les variations de chaleur et de temperature dQ et dT sont liees par 
la relation suivante 

dQ = c dT 

ou c est une constante appelee chaleur specifique de la transformation poly- 
tropique. On voit que les transformations adiabatiques sont des cas parti- 
culiers de transformations polytropiques avec c = 0. 
Les gaz parfaits obeissent a l'equation d'etat (Jlj) 

T 

P(q) = —pm{q) 

m 

qui relie la pression P, la temperature T, la densite de masse p m et la masse 
m d'une particule du gaz. Pour un gaz parfait subissant une transformation 
polytropique, la relation entre la pression P et la densite de masse p m est 
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de la forme [Sj 



p = k p j (A.r 

ou K est une constante et ou le coefficient polytropique 7 est 



Cp c 



7 = , (A.2) 



c v c 



Cp et c v etant respectivement la chaleur specifique a pression constante et la 
chaleur specifique a volume constant. Pour les transformations adiabatiques 
(c = 0), on retrouve la valeur bien connue du coefficient adiabatique 7 = 

c iM- 

En outre, en utilisant l'equation des gaz parfaits (DO) et en introduisant 
la temperature polytropique T 7 , temperature pour laquelle la densite de 
masse p m est egale a 1, on trouve que la constante K est egale a 

T 
m 

Nous allons voir maintenant que les etoiles sont des gaz autogravitants 
polytropiques. 



A.2 Les etoiles 

Les reactions thermonucleaires qui se deroulent au coeur des etoiles con- 
stituent leur source de chaleur. Celle-ci se propage du centre chaud vers 
la peripheric plus froide puis est rayonnee a l'exterieur des etoiles. II est 
raisonnable de supposer que les etoiles sont en equilibre convectif, ce qui 
veut dire que les transferts de chaleur du centre chaud de l'etoile vers la 
peripheric plus froide se font par convection, les transferts par conduc- 
tion etant negligeables devant ceux-ci. Ainsi chaque element de gaz con- 
serve sa chaleur et se transforme done adiabatiquement dans l'etoile. Ainsi 
les etoiles peuvent etre considerees comme des gaz autogravitants poly- 
tropiques. Nous allons maintenant determiner l'equation de la densite des 
gaz autogravitants polytropiques en equilibre hydrostatique. 
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A. 3 Equilibre hydrostatique 



Nous allons determiner l'equation de la densite d'un gaz autogravitant 
polytropique en equilibre hydrostatique et obeissant a l'equation d'etat 
(lA.ip . En utilisant cette equation et l'equation ( 11.51 ). on obtient 

Vq(Pm ) = Tf, Pm{q) ■ 

En introduisant l'indice polytropique 

n = ^- = ^— -, (A.3) 

7-1 c p - c v 



en posant pour la densite 



Pm = Po ^ 



ou p est une constante, et en introduisant le rayon vecteur sans dimension 
A defini par 



q = a A , a 



T 7 (n + l)p " 1 



4tt G m 

on trouve l'equation des systemes autogravitants polytropiques 

v\e + e n (\) = o . (A.4) 

Dans le cas de la symetrie spherique, cette equation devient 

hU*£) + e '=°- (A - 5) 

Cette equation qui est l'equation de la sphere polytropique est appelee 
equation de Lane-Emden d'indice n. Si l'indice est le meme dans toute 
la sphere, on peut poser que p est la densite au centre et en deduire la 
premiere condition initiale a savoir 

0(A = 0) = 1 . 

Pour que l'equation soit reguliere en 0, on impose la deuxieme condition 
initiale 
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En la resolvant, on en deduit la densite et toutes les grandeurs physiques 
comme la temperature et la pression, comme pour la sphere isotherme. 
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Annexe B 
Theorie de Jeans 



La theorie de Jeans [21 [2BI EH] montre comment se forment des conden- 
sations de matiere a partir d'un fond de matiere homogene. Cette theorie 
dont le but est d'expliquer la formation des galaxies est particulierement 
interessante par sa simplicity. Lorsque de petites variations de densite se 
propagent sinusoi'dalement, le systeme est stable. Par contre, lorsque ces 
petites variations croissent exponentiellement, le systeme devient instable. 
Determinons tout d'abord l'equation de propagation de ces petites pertur- 
bations. 



B.l Equation de propagation 

Soit un fluide autogravitant de densite de masse p m , de pression P, de 
vitesse v et creant un champ de gravitation g. Le fluide est regi par les 
equations de la mecanique des fluides (equation de continuite et equation 
d'Euler) 



%1 + V.(^)=0 , f+^x v)v = -±-VP + S (B.1) 
et les equations de la gravitation newtonienne 

Vx^=0 , V .g = -4nG p m . (B.2) 

Les perturbations au premier ordre sont determinees par rapport a un 
fluide statique uniforme ou les effets de la gravitation sont ignores. Pour le 
fluide statique, on a 
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Pm — Po = constante , P = Po = constante , u = , g = . 

Considerons une perturbation de ce fhiide statique uniforme. Soit respec- 
tivement pi, Pi, v\ et g\ la densite de masse, la pression, la vitesse et le 
champ de gravitation de cette perturbation. Au premier ordre, les equations 
de la mecanique des fluides et les equations de la gravitation newtonienne 
deviennent 

^ + Po V.ff 1 = , § = -lvF 1 + ?1 
at ot po 

et 

Vx^ = , V.g=-4<irGpi. 
On introduit la vitesse du son du fluide v s [2H] 

dP Pi 



v s = 



9pm Pi 

Combinant ces equations, on trouve l'equation de propagation suivante 
pour la densite 



d 2 P i 



v 2 s V 2 pi + 4 7T G po pi . (B.3) 



dt 2 

Nous allons en deduire la relation de dispersion de ces petites perturbations. 

B.2 Relation de dispersion 

Les ondes planes 

pi(q , t) a exp [i(k .q — to t)] 

sont solutions de l'equation de propagation avec la relation de dispersion 
entre la pulsation u et le vecteur d'onde k 

u 2 = k 2 v] - 4 7T G po . (B.4) 
On introduit le vecteur d'onde de Jeans 
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V4 7T G po 
kj = 

v s 

qui est le vecteur d'onde pour lequel la pulsation u ( IB. 41 ) s'annule. En 
introduisant la longueur d'onde A = -M 1 et la longueur d'onde de Jeans 



on a 



co 2 = 4 tt v 2 ^_ - _ j . (B.6) 

Pour les longueurs d'onde plus petites que la longueur d'onde de Jeans 
(A < Aj), on a u 2 > 0. La perturbation varie sinusoi'dalement. II n'y a 
pas formation de condensation. 

Pour les longueurs d'onde plus grandes que la longueur d'onde de Jeans 
(A > Aj), on a u 2 < 0. La perturbation croit exponentiellement avec le 
temps. II y a done formation de condensations. 

Nous allons maintenant voir le cas des fluides homogenes isothermes. 



B.3 Instability dans les fluides homogenes isothermes 

Des condensations se forment seulement a partir de perturbations ayant 
atteint une taille critique qui est la longueur de Jeans Aj. Des perturbations 
spheriques se condensent si leur rayon Q est superieur a 4f . Considerons 
un milieu homogene isotherme compose par un gaz parfait de particules de 
masse m a la temperature T. La vitesse du son au carre est v 2 = — et la 

masse de la perturbation spherique de rayon Q est M = m N = pq. 
D'apres la relation ( IB. 51) . on en deduit, dans ce cas que la valeur de notre 
parametre r) R = G q T n ( 11.171) est 



4 7r G po Q 2 4 (tiQ 



vl 3 V A 



Des condensations spheriques se forment lorsque son rayon Q est superieur 
ou egal a la moitie de la longueur de Jeans Aj. D'apres la relation ( 1B.7I ). 
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la valeur de 77 a partir de laquelle la perturbation se condense est 



Pour 77 < 77 j, la perturbation oscille sinuso'idalement, elle ne se condense 
pas. Pour 77^ > 77J, la perturbation s'effondre sur elle meme , il y a 
formation de condensation. La theorie lineaire de Jeans donne pour 77^ 
la valeur d'instabilite r]j (IB. 8ft . La valeur du parametre 77^ pour laquelle 
la sphere isotherme collapse dans l'ensemble canonique est 77 ca „ = 2.43... 





(eq.(103)). 
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Abstract 



The self-gravitating systems are formed by particles interacting through gravity. They 
describe structure formation in the universe. As a consequence of the long range inter- 
action of gravity, they are inhomogeneous even at thermal equilibrium. They can be in 
a gaseous phase or in a collapsed phase. We formulate the statistical mechanics of the 
self-gravitating systems. The thermodynamic limit where the number of particles N and 
the volume V tends to infinity with fixed is relevant for the gaseous phase. The do- 
mains of stability of the gaseous phase are different in the microcanonical ensemble and 
in the canonical ensemble. The instability of the gaseous phase leads to its collapse into 
a phase of infinite density. In the thermodynamic limit, the mean field approach gives an 
exact description of the gaseous phase. After introducing the self-gravitating systems with 
one kind of particles (chapter 1 and 2), we study the self-gravitating systems with sev- 
eral kinds of particles (chapter 3) and the self-gravitating systems in the presence of the 
cosmological constant (chapter 4). We formulate for these two types of self-gravitating 
systems the statistical mechanics and the mean field approach describing the gaseous 
phase. We find the equation governing the density of particles. We explicitely compute 
thermodynamic quantities and find that they are extensive (proportional to N). We ob- 
tain the domain of stability of the gaseous phase. In the self-gravitating systems with 
several kinds of particles the density of the light particles is flatter than the density of 
the heavy particles. Scaling exponent of the self-gravitating systems with several kinds of 
particles are computed. The cosmological constant acts as an uniform density of energy 
with a repulsive gravitational effect on the matter. The particle density is a decreasing 
(increasing) function of the radial distance when the self-gravity dominates over the cos- 
mological constant (and vice- versa). Monte Carlo simulations show that the mean field 
describes the gaseous phase with an excellent accuracy. They allow to study the collapsed 
phase and confirm that the phase transition happens when the isothermal compressibility 
diverges. The presence of the cosmological constant extends the domain of stability of the 
gaseous phase. 
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Resume 



Les systemes autogravitants sont constitues de particules interagissant mutuellement 
par la gravite ; ils decrivent la formation de structures dans l'univers. Comme consequence 
de l'interaction a longue portee, les systemes autogravitants ne sont pas homogenes meme 
a l'equilibre thermodynamique. Ils peuvent exister sous une phase gazeuse ou sous une 
phase collapsee. La mecanique statistique des systemes autogravitants est presentee. La 
limite thermodynamique ou le nombre de particules N et le volume V tendent vers l'infini 
avec fini est pertinente pour decrire la phase gazeuse. Les domaines de stabilite de la 

phase gazeuse sont differents dans l'ensemble microcanonique et dans l'ensemble canon- 
ique ; l'instabilite de la phase gazeuse entraine son collapse dans une phase de densite in- 
finie. L'approche du champ moyen de la mecanique statistique decrit exactement la phase 
gazeuse dans la limite thermodynamique. Apres avoir presente les systemes autogravitants 
ne comportant que des particules identiques (chapitre 1 et chapitre 2), nous avons etudie 
les systemes autogravitants comportant plusieurs sortes de particules (chapitre 3) et les 
systemes autogravitants en presence de la constante cosmologique (chapitre 4). Pour ces 
deux types de systemes autogravitants, nous avons developpe la mecanique statistique puis 
nous avons developpe l'approche du champ moyen decrivant la phase gazeuse. Nous avons 
trouve l'equation verifiee par la densite de particules, nous avons explicitement calcule 
les grandeurs thermodynamiques et nous avons montre qu'elles sont extensives (elles sont 
proportionnelles au nombre de particules N). Nous avons determine le domaine de sta- 
bilite de la phase gazeuse. Dans les systemes autogravitants comportant plusieurs sortes 
de particules, la densite des particules legeres est moins contrasted que la densite des 
particules lourdes. Nous avons calcule les exposants d'echelle des systemes autogravitants 
comportant plusieurs sortes de particules. La constante cosmologique agit comme une 
densite d'energie uniforme ayant un effet gravitationnel repulsif sur la matiere. La densite 
de particules est une fonction decroissante (croissante) de la distance radiale lorsque l'au- 
togravite domine la constante cosmologique (et vice- versa). Nous avons effectue des calculs 
Monte Carlo pour les systemes autogravitants en presence de la constante cosmologique. 
Ils permettent d'etudier la phase collapsee ; ils confirment que la phase gazeuse est decrite 
avec une grande precision par le champ moyen et ils montrent que la transition vers la 
phase collapsee s'opere lorsque la compressibilite isotherme diverge. La presence de la 
constante cosmologique etend la stabilite de la phase gazeuse. 



